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BTROrUoTIOlN 

This  paper,   an  expansion  and  discuesion  of  a  series  of  lectures  given 
by  Sr.    'f.  K.    bcheubel  cf  the  Technical  Institute,    Darraetadt,   Germany,    was 
undertaken  by  the  writer  as  a  thesis  project  while  working  for  the  degree 
of  Master  of  Science   in  Aeroneutical  Engineering  at  the  University  of  Kich- 
igan,   Ann  Arbor,   Michigan.      These  lectures  were  given  at  the  University  of 
Michigan  during  the  fall  of   1951  while  Dr.   Scheubel  was  visiting  the  Uni- 
versity at  the   invit£.tior   of  Dr.   F.   W.    Jonlon   of  the  Aeronautical  Engineer- 
ing'Departaent. 

Dr.   Scheubel  presented  the  material  in  six   lectures,    the  first  three 
devoted  to  quasi-static   stabilitv  concopte,    the    last  three  to  their  use   in 
the  solution  of  the  equations  of  motion.      His   limited  time   prevented  a  de- 
tailed accovmting  of  the  assunptions  involved  as  well  as  minute  explanations 
of  the  approximations  made.      It  is  the  purpose  here  to  verify  and  expand  ,on 
his  presentation  to  a  degree  consistent  with  the   limited   scope  of  such  a 
paper. 

Briefly,    his  material  covered  the  folloving.      His   lectures  were  re- 
stricted to  symmetric  or  longitudinal  motion.     He  developed  the   longitudi- 
nal equations  of  motion  and  the  quasi-etatic   stability  criteria  at  equili- 
brium and  at  constant  speed.      He  solved  the  equations  using  these  quasi- 
static  stability  concepts  for  both  the  phugoid  and  short  period  modes  for 
the  stick-fixed  case.     An  anplltude-phace  relation  for  the  two  variables 
applicable  to  the  Dodee  was  discussed.      He  then  introduced  the  degree  of 
freedom  about  the  elevator  hinge  line  and  solved  for  the  stick-free  case. 
Finallv  the  effects  of  an  elevator  impulse  were  discussed. 

This  approach  to  dynamic    longituuinal  motion  differs   fx^om  standard 
methods  mainly  in  the  quasi-static   stability  concepts  as  regards  their  in- 
sertion in  the  solutions  to  the  ecuations  of  motion.      In  this  respect,   the 
method  may,    as   illustrated,   be  safely  applied  only  to  conventional  air- 
craft,   i.e.,   a  rigid  body,   where  the  number  of  degrees  cf  freedom  and  thus 
the  complexity  of  solution   is  restricted,      rhe  analysis  of  the  motion  as- 
Eunlng  a  phygold  and  a   short   period  oscillation   serves   only  to  illustrate 
the  handling  of  the  quasi -static  concepts.     The  limitation  of  the  method 
in  respect  to  acceptable  separation  of  the  motion  into  these  modes  was 
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broxight  out  by  Dr.    Scheubel  and    is  discuEBed.      The  form  of  the   eolution   of 
the  complete   set  foregoing  this   eeparatlon   is   included   in  this   paper. 

The  equations  of  motion  as  eetfiblished  "by  the  standard  approach  of 
reference  1  are  relegated  to  the  flight  path  direction  or  wind  axes  in 
Older  that  the  further  discussion  vill  bo  based  on  common  aeBumptione 
and  approximationB.  Dr.  Scheubel 's  equations  as  developed  for  a  eyetera 
reletivc  to  the  flight  path  v;ill  then  be  ccr.pared  term  by  terra  with  the 
standard  equations.  Frora  this  point  onwsrd,  i.e.,  beginning  vith  non- 
diraensionalizing  the  equations,  the  presentation  is  as  Dr.  Scheubel  de- 
veloT)ed   it. 


koKenciature 


m  =  W/g 
u',   w 


X,    Z,   M 

e,  (f-  ,  i?,    or  ££,  L^  ^  ,h  y< 


I       =   mi 

yy-        y 

oil 


V  =  V      +  AV 
o 


L,    D,   M 
W 


o'      o 


'L 
P 


S 


u  =  AV/V 


mass  of  the  aircraft 

initial  velocities  along  i  and   ^  body 
axis 

perturbation  velocities  along  x  and   z 
body  axis 

angular  velocity  about  y  axis  (body  or 
vind) 

air  forces  and  moments  * 

perturbation  of  attitude  angle,    flight 
path  angle,    and  angle   of  attack 

initial  attitude  angle 

aircraft  radius  of  gyration  about  y 
axis 

moment  of  inertia  about  y  axis 

initial  flight  path  velocity 

velocity  pert lurbat ion  along  flight  path 

velocity  along  flight  path  folloving  per- 
turbation 

lift,   drag  and  air  moment 

aircraft  weight 

initial  thrust  and  drag 

elevator  deflection 

lift  coefficient  for  aircraft 

air  density 

wing  area 

aircraft  drag  coefficient 

non-dimensional  velocity  perturbation, 
ratio  of  perturbation  tc  total  flight 
path  velocity 


e   g  ^  VS 


M   = 


2W 


6  =  t/t„ 


M 


EX  il 


y    y 


*r*2 


n 


'Lh 


h 
Go 


T  ,  T 
p'   r 


P 


or  ApAg 


adopted  standard  time  unit 

mass  density  (See  Appendix  A) 

non -dimensional  time 

initial  moment  about  y  axis 

tail  length,  distance  from  aircraft 
e.g.  to  center  of  pressure  of  horizon- 
tal tail 

con.^ugate  complex  root  of  characteristic 
equation.  E  is  real  part,  I  is  imaginary 
part 


constant  coefficients  of  equations 

position  of  aircraft  center  of  gravity 
from  most  forward  part  of  aircraft  in 
percent  of  wing  chords 

position  of  neutral  point  of  aircraft 
from  most  forward  part  of  aircraft  in, 
percent  of  wind  chords 

angle  of  attack  of  hori. ontal  tail 

lift  coefficient  of  horizontal  tail 

area  of  horizontal  tail 

angular  frequency  of  mode  of  motion 

periods  of  phygoid  (p)  and  rotary  (r) 
modes  of  motion 

wave  length  of  phygoid 

distance  "between  the  hinge  line  and  cen- 
ter of  gravity  of  the  elevator 

mass  of  the  elevator 

radius  of  gyration  of  the  elevator 

hinge  moment  coefficient  of  elevator 

mean  aerodynamic  chord  of  the  horizon- 
tal tail 


Ir  .  elevator  impulee 

H  constant  altitude  reference 

o 

h  altitude  change 


Subscript  «  indicates  eteady  state. values  in  the  response  discussion. 
Any  other  symbols  used  are  either  obvious  or  are  locally  defined  for  ease 
in  handling  equations. 


II.   Development  of  the  Method. 

A.   The  eqxiations  of  longitudinal  motion. 

The  standard  development  in  a  perfectly  general  vay  es  accompllBh- 
ed  in  reference  1  result  in  the  equations  of  longitudinal  motion  of  the  form, 

A.l   m(u'  +  W.q)  =  X  ,u'  +  X  v  +  X  q  -  mg  cos  ♦en,©   .  I 
1^     u      w     q^    ^       1    / 

A. 2   m(v  -  U-q)  =  Z  ,u'  +  Z  w  +  ?  q  -  mg  ein»€H,9 
1      u'      w     q^   ^       1 

A. 3       mi~q=Iq=Mu'+Mv+Mq 

y    yy    u     w    q 

These  equations  are  relative  to  body  axes  and  contain  the  follov- 
ing  aesumptlons  and  approximations: 

(a)  Initial  symmetric  steady  motion  is  assumed. 

(h)  The  air  reactions  do  not  depend  on  the  rates  of  change  of  the 
variables,  U,  V  and  W,  or  their  integrals. 

(c)  Second  order  and  higher  terms  of  the  air  reactions  are  neglected, 
i.e.,  only  infinitesimal  disturbances  from  an  initial  steady  motion  are 
treated. 

(d)  The  aircraft  has  a  plane  of  symmetey  and  the  steady  motion  about 
vhich  the  disturbances  occur  is  symmetrical  with  regard  to  that  plane. 

(e)  The  disturbance  initially  imposed  on  the  system  is  unenforced  and 
the  controls  are  locked.  This  rules  out  an  initial  couple,  M  ,  and  the  vari- 
ation of  M  vlth  S   . 

The  variables  are  seen  to  be  u,  w  and  q.  These  are  illustrated  in  Figxire 

1.  The  relation  between  the  body  axes  and  the  wind  axes  is  shown  in  Figure 

2.  The  body  axes  are  primed.  Angles  are  measured  positive  counterclockwise 
from  the  horizontal  reference  for  T  and  ^^^-i,  from  the  wind  line  f or  o(  . 

In  regard  to  these  figures  it  is  seen  that, 

(a)   Since  q  is  a  velocity  and  9  a  displacement,  6   =  /q  dt  =  A  C»(  -i-  A  ^ 


ijofiiZOsTAL 


(b)  V  is  the  steady  Etate  velocity  and,  V  =  V/t  -i-  Ut. 
Then  there  Is  the  definition,  V  =  V  +  AV 

'     Gt    dt 

(c)  In  essence  u  is  nov  iJ.V,  and  w  does  not  exist. 

The  inei^la  terns  of  the  equations  of  notion  parallel  and  perpendicular 
to  the  wind  direction  and  aljout  the  center  of  gravity  become,  neglecting  the 
Inertia  force  due  to  linear  acceleration  perpendicular  to  the  wind  direction, 


for  A.l: 
for  A. 2: 
for  A.3: 


m 


d(/^V) 


dt 


dt 


d^(  o<  ^  y  ) 


77 


dt" 


The  weight  components  relative  to  the  wind  axes  "become, 
along  x:    mg  sin  ^ 
along  ^;   -mg  cob  dy 
and  these  are  only  affected  by  A  l^  . 

There  remains  only  the  air  reaction  derivatives.  The  force  created 
on  the  tail  due  to  q  is  the  largest  resulting  from  this  distur-bance  and 
from  experience  it  is  known  to  be  small  along  the  wind  axes  so  X  and  Z 

Si        H 

are  neglected. 
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In  the  standard  method  the  left  eutid  drag  contributed  to  the  air  reac- 
tion forces  "both  along  the  chord  and  normal  to  it.  Also  velocity  pertur- 
bations existed  along  both  axes.  L'sing  vind  axes,  only  the  velocity  per- 
turbation, AV,  nov  exists.  Any  sinking  velocity  along  the  z  vind  axis  is 
merely  a  change  in  angle  of  attack.  The  air  reactions  along  i  consist 
only  of  the  drag  in  its  entirety,  neglecting  variations  of  thrust.  The 
drag  varies  with  both  <^  and  V.  The  air  reactions  along  z  cOTisist  cnly  ' 
of  the  lift  in  it£  entirety  and  is  affected  by  the  sane  quantities.   A 
pure  rotary  perturbation  about  the  center  of  gravity  changes  both  o<  and 

(J^but  only  CX  affects  the  linear  forces.  The  same  applies  to  the  air  nio- 
rnents.  The  variations  with  <^  are  found  by  vind  tunnel  tests  and  any  rota- 
tion of  the  r.iodel  about  the  center  of  gravity  ie  considered  pure  ^    .      This 
does  not  hold  for  the  case  of  the  effect  on  M  of  both  ^  c<  and  l^  ^  .     These 
two  perturbations  make  up  6,   which,  when  multiplied  by  the  tail  length, 
gives  the  effective  sinking  velocity  of  the  hori^ontal  tail.  This  sinking 
velocity  in  turn  is  felt  as  a  change  in  angle  of  attack  of  the  tail. 
Thus  the  external  forces  for  A.l  become, 

^  ^91  .  5^  ^7  -.  ^-^  (mg  em  J-  )  AT 

and  equation  A.l  is. 


d(^V)       SB      _        Sd  ,„  v<  A  v^ 

"^       dT  =  S^  -<^    -  §v  '^V  H-   tng  cos  J^  A  r 


Equation  A. 2   is, 

"^"^"^~%l^°^    "•"SAT'-^V+mg  sin  2f-  A  (T     .,. 
Equation  A. 3  becomes, 

T  d-^(  o<.       (T  )  ,  SM  ,_,  hV.   ,,,  5K   .V  bM      ^ 

I  * r *-  =  r-r-  AO*'  ^   ^rr  ^Y  +   ^-j-  _o<  r    v-r  ^^  dr    . 

i                     w            ,.2  o^  OV  o«  or- 

■  " "  at 

The  use  of  the  wind  axes  then  has  the  advantage  of  the  aircraft's  for- 
ward motion  being  along  the  i  axis  so  that  the  lift  and  drag  forces  always 


lie  along  these  axes  and  reeoluticjn  into  nonnal  end  chord  forces  ie  not 
neceesary.  Furthermore,  no  resolution  of  the  velocity  is  required  to  be 
made  along  the  ajcie  of  perpendicularity  vhich  le  the  case  if  the  forces 
are  resolved  along  the  body  axes.  The  disadvantage  ie  that  these  axes 
change  continuously  eo  that  the  moments  and  products  of  inertia  change. 
For.  BBiall  disturbp-nces  from  initial  hori<:ontal  flight  these  diecropancieB 
are  considered  negligible.  Also  the  thrust  does  not  act  along  the  vind 
axes  for  all  aircraft  and,  for  an  individual  aircraft,  at  all  times. 

The  continuity  of  this  payer  and  the  advantage  of  physical  "feel" 
that  is  had  from  the  development  of  the  equations  of  motion  from  initial 
suppoeiti'on  of  wind  axes  can  be  beet  maintained  by  including  such  a  de- 
velopment in  its  entirety  at  this  point.  This  is  the  method  used  by  Dr. 
Scheubel  and  appears  to  be  in  general  use  in  Germany.  See  Eeference  2. 


Figure  3. 


With  the  aircraft  in  an  attitude  and  with  the  external  forces  as 
shown  in  Fit^ure  3>  the  external  forces  are  equated  to  the  time  rate  of  change 
of  momentum  In  the  direction  of  the  wind  axis,  i.e.,  along  V,  and  one  equa- 
tion of  equilibrium  follows. 


^     ^  =  T  -  D  -  W  sin  ^'   . 
S  dt 


IC 

The  external  forces  are  ceueed  to  chajige  by  8m£ill  changee  or  pertur- 
bations in  the  valueG  of  cK    ,  V,    ^T*  ej\±    T  .  The  forces  in  equillbritm 
then,  follovlng  a  small  disturbance,  are  the  initial  values  plus  the  in- 
crements  due  to  the  changes,  \ 

^  ^  =  T  -  r  -  V  sin  Jf   -  II  ^V  -  II  ^V  -  1^  .^cx  -  W  cos  Jr  A  ^  .  \; 
S  dt    o    o  o   oV     oV     dJ^  V 

This  is  the  result  of  a  Taylor  Series  expcjneion  of  the  quantities  T,  5 
and  W  as  functior.s  of  «^  ,  V,  ^  and  ^   in  vhich  second  order  terms  and  a- 
bove  are  neglected  under  the  assumption  that  the  changes  are  small.  Assuming 
an  initial  condition  of  equilibria^,  it  is  seen  that  this  eqijation  is  compar- 
able to  equation  fi.l   vith  one  exception.  Thus  all  the  assumptions  previously 
enumerated  h<?ld.  The  exception  is  the  thrust  term.   It  is  easily  seen  that 
thrust  is  not  affected  by  chajiges  in  0<  ,  iT*  ,  or  ^  .  The  stick-fixed  case 
is  considered  here  also.  Thus  all  terms  involving  ^  are  neglected. 

This  first  equation  is  nov  manipulated  into  a  form  necessary  for  sub-  . 

Sequent  solution.   Initial  equilibrium  and  assvuning  JT'  smell  gives, 

T  -  D  -  V  sin  r   =  0   and,    cos  T  =  1   or   W  =  L 
o    o  o  ' 

also,   L  =  C,  ^  fv^S   and,   D  =  C^^  ?   V^S 

*"        T  L 

o     V  o     V  there  is 

Multiplying    terms  in  T  by  =f-  and  -  ,  and  terms  in  D  by  —  and  - 

o  o 

^  ^  .v    I-    ^'^    tl      r  1-  ^^  t^   ^^  -/v   T   ^ 
g  dt  "  o  T   5v   V  "  o  D   5?   V  "  o^  ^"  ' 
o  o 

^'■°^'        _v'   ^T       ^    In  T    .   _V     hD       S    In  D       ^""^   letting       ^ 

T  .    §V  "  d    In  V    '   D       Sv       ^    In  V  V  '    ^ 

o  o 

there  is, 

w    dv_  V    ^^Vl^_  w    d(.:.v) 

g     dt       g  dt  "   g       dt 


2 


Dein{3  "the   standard  notation,  ^~rr  =   --  and  dividim^  bv     P  :r-  S,  there   is 
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5  p;  ir^=  --«  -«  -  - -^  -  -  ^s-iTv  -  r  rv7^  -  • 

Thie  eqiiation  Is  non-dineneiane  li„ed   in  the  ueual  ia3.r.ner.      The  right 
side  of  the  equation  is  non-dimensional  as   it  stands.      Considering  the   left 
side,   multiplying  nunerator  by  V/V  and  the  denominator  "by  t  /t  , 

W    2    V_  d(a'/v) 


3 


P  v's 


T 


W    "^    V 

A  tine  unit,  t  ,  is  adopted  "bv  lettin(^,  ~  —    —  =  1  then, 

'^  S  ^y-S  tg 

1/^       l/^ 

t  =  — ■= xTl  -  M  -Tj —  =  C.  ~  vhere  the  mass  nu.iiber,  M    ,    is  given  "by, 

^       g^V^  /-    /        V  L  g  / 

yti   =      -  ■   '/■    =  (~)  ■\j\   .      This  mass  number,    or  density  factor,    is 

discussed   in  Appendix  A.     I'ov  the  non-dimensional  time   is  given  by,    t/t     =  C 


g 


60  that,  d{— -)  =  d  d.  and  the  variable  lN   is  transformed  as,  —•  =  u. 

t  / 

g 

The  first  equation  of  motion  then  becomes, 

A   ).  du        .  ^  «/         -^         >*       -     /^    IJi  r.       "^o    SlnT. 

d?^  L^  ov  o    o    In  V       i>       o    Jji  V 

*^  o 

Again  referring  to  Fi^iure  3;  the  external  forces  perpendicular  to  the 
vind  line  are  equated  to  the  centrifugal  force  or  mas^  times  the  centrifugal 

acceleration.  -  N  -rr  -   -  Y  rrr  -   L-W  cfia  r    since,  -  \  ^-r  =  C,  because 

g   at   g   at  '  3   at 

do( 
V  77  results  in  no  change  in  velocity  perpendicular  to  the  flight  path  and  the 

linear  acceleration  in  this  direction  is  neglected.   If  a  small  disturbance 

Is  introduced,  and,  since  lift  does  not  vary  with  ^    ?.nd  the  etick-ifixed  case 

is  assumed. 


-  V  r7'^=    L      -  W  cos    r      ^  ?3  -C<     '   ?T7  -V  -  W  sin  J^  A  T  .  a 

~       dt  o  o      OiX  dV  , 


With  the   same  assumption  as  before,    that   of   initial  equilibrium,    this 
equation  is   seen  to  be  the  seme  as  equ-Ttlon  A. 2  and   so  the  same  assumptions 
and  approximfltlor.s  are   implied.     Expanding  the  derivatives, 
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and 


nr> 


EO 

2  a    ]J1  C 

/aL  =  P  ~  i3  (:,        ^(X   +   G^    (2    ^     ^    ,,  -ir)    u   ) 


P^^c^^^^  .0^(2  ^  -y-3^; 


alto,    W  Ein  ^  ^  C"  =  C-    sir    i^   P   —  S  ._i  S**  ,    so  th-it  the  entire  equation  "be- 

CGmes, 

y         ^  ,-.  ^    In  C^ 

"     pv  S       ''^  ^CX  -  d    Ir.  V  L 

The   left  side   ie  r.on-dimenBion&lized  ae  "befoie.     Cividing  the  demommr  tor 

'/ 

by  tVt^,    t'^-cre   is, 

.     s  , 

V     2       v_     1  y*   ^  d_J-  ^    • 
^pr^s    ^s   d(|-)     ^'^ ., 

So  finally  there  is  the  second   equation  of  notion, 

:  .  a    In  J^ 

A. 5  ^  -  J^    (2    ;   s-r-rr^)   u  -  C^       ^(X    -  J,    Gin  ^  .^  iT  =  0   . 

The  final  equation  of  the  longitudinal  motion  is  oht^^-ined  "by  equating  the 
motnentt.  to  the  moment  of  moment  urn  about  tbe  center  of  gravity  of  the  aircraft, 
oonsidering  Figure  3  again,   thie   ie, 

W  ,2  d-(  (X       !^)       ^^ 

^  y-    dt^  ° 

W  2  2 

where,  -  i  is  the  moment  of  inertia  rnd  1  is  the  radius  of  .nri^ation.   Since 
6  7  y  ^^ 

M  i£  a  functiori  of  lift  and  drag  which  in  turn  are  not  affected  by  tf*  ,  nejlect- 
ino  the  influence  of  V,  .^nd  vith  the  stick-fixed  nttum.  tior ,  there  ie, 

E  i-  d"(  Q^   ^  )  ^  J       d^(  g^  ir  ]  ' 

8  '^      dt-    ^  3ry     ^t  . 

=  ^o  ^55^-^  ^  5v-^  -  5^  ^   r^-^- 

This   Ie   seen   to  be   the   ecme  ae  equcition  ."• .  3  aseuralns   initial  equilibrium. 
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Putting  the  change  in  nonent  in  coefficient  form, 

u       ^^M  o  V^  e3/2  .,^   V  "^^W  O  V^  ^3/2  ^V  ^  ^-^M   p  V^  „3/2  d0 

^^dt^ 

Since,. c^.Af=^i^,^=||. 

An  important  point  to  note  here  is  the  Gemu-n  uee  of  the  square  root  of 
the  wing  erea  ae  the  additional  length  term  when  converting  moments  to  coef- 
ficient form.  The  reasoning  is  that  this  quantity  is  more  easily  definable 
for  the  variegated  types  of  wings  now  in  existence  whereas  the  m&n  aerodynamic 

chord  is  somewhat  nebulous  in  definition  in  the  literature.   Another  reason 

1/2 
le  the  fact  that  the  qi;iantity  rj^^'~   ranges  from  C.8  to  1  IC  end  this  facili- 
tates its  approximation  to  unity  where  it  appears. 

Considering  the  last  term  and  non-dimensionfiliiing  the  derivative  bv  lett- 

da  .^^V2 

ing,  q  =  TT 


dt   V 


^^M      p  a  ,3/2  d0       ^^M  p  £  „3/2  de      l}lz. 
Then  non-dimensionalizing  the  entire  term, 

Now  nan-dimfcnBionaliz.ing  the   left  side  of  the  equation,   after  dividing  bv 


Wl2 

g      V 


2  1      i^(o<      r  )      ,    ,2    iv-_.;j      li. 


.2 

This  equation  nay  be  rewritten  in  the  more  concise  form, 

..tic-  C„  §^CX-C    -|/,-.«    ^.,^     i^i-^-^     ^^u-O 
CI        "^  i  q  1  i 

y         7  7  7 

There  are  then,    the   longitudinal  equationc  of  motion  in  nor-dlmeneicn&l  fort 


lu 


3  In  D   ^o  ^  In  T. 


^■^        "^^D'rrrf-r  s-^'^-'^r.-"    x-^'" 


o 
^  In  G, 


Ql  0\     1     '  01  1 

y        y  y  y 

These  equations  have  been  seen  to  comoare  term  by  tern,  except  for  the 
J.nclusion  of  variation  of  thru£t  with  velocity,  with  theee  developed  in  a 
more  general  sense.  The  aeeumptions  and  approximatione  have  been  cited.  The 
dependent  variables  are  u,  C-  and  /^ . 

A   more  general  case  would  be  that  of  aesuming,  instead  of  the  homogeneous 
set  above,  the  equating  of  the  quantities  above  to  forcing  functions  of  the 
elevator  displacement,  f  .     This  would  require,  in  addition,  en  equr^^tion  of 
motion  involving  moments  about  the  elevator  hinge  line.  This  would  be  the 
case  of  stick -free  stability  ajid  motion. 

It  is  noted  that  a  striking  difference  between  these  equations  and  those 
more  fam.liar  to  most  students  is  the  form  in  which  the  stability  derivatives 
have  resulted.  Since  the  problem  is  now  ready  for  analysis  of  modes  of  motion, 
damping  factors  and  associated  desired  results,  the  quest  on  arises  as  to  why 
this  form  of  the  derivatives  and  how  they  are  to  be  evaluated.      ^ 

The  usxial  procedure  is  the  separate  evaluation  of  each  derivative  by 
similarity  to  previous  determinations  or  by  wind  tunnel  or  flight  testing.   It 
is  seen  that  certain  of  the  derivatives  in  the  three  equations  above  are  in 
combination.  For  a  special  but  usual  case  of  the  solution  of  the  equations  of 
motion  further  combinations  of  these  derivatives  appear.  A  discussion  of  both 
the  general  case  and  this  special  one  of  the  solution  foliows. 
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B.   Solution  of  the  Equations  of  Motion 

The  most  general  approach  to  the  solution  of  the  equations  of  motion  ■be- 
comes quite  involved  in  that  the  solution  is  carried  to  the  determination  of 
an  equation  for  each  variable  composed  of  both  the  complimentary  function  and 
particular  integral.  See  Beference  1.   If  the  eqxiations  are  reetrictid  to  the 
stick-fixed  case  and  the  motion  is  atsumed  resulting  in  the  ab-eence  of  an 
applied  forcing  function,  the  equations  will  be  homogeneous  and  the  solution 
resulting  will  be  the  complimentary  fvmction  only. 

The  discussion  here  will  be  further  restricted  in  that  the  equations  in 
terns  of  the  variables  will  not  be  the  end  result   This  would  merely  give 
the  particular  motion  resulting  from  certain  initial  conditions. 

The  characteristic  equation  which  results  from  the  solution  of  the  dif- 

7  Z 

ferential  equations  being  of  the  fona,  x  s  x  e  ,  will  be  analyzed  for  Ite 
roots  and  information  relative  to  inspection  of  those  roots.  This  will  show 
the  character  of  the  motion  as  regards  stability,  periodicity  and  demplng. 

The  nature  of  the  characteristic  equation,  or  stability  quart ic,  wij.1  be 

set  down  from  the  determinant  resulting  from  the  solution,  x  =  x  e      Tx- 

o 

perionce  has  shown  that  the  two  sets  of  roots  that  form  the  solution  are 
usually  of  such  widely  separate  magnitude  that  they  may  be  separated.  The 
special  case  mentioned  is  based  on  this  fact 

Since  Scheubel's  solutions  employing  the  quasi -static  stability  concepts 
are  based  on  the  premise  that  the  roots  are  separable,  the  general  characteris- 
tic equation  will  be  set  down  without  further  comment  and  the  special  case 
will  bo  stressed  in  preparation  for  the  solution  by  his  quasi-static  stability 
concepts.  The  errors  involved  in  this  approximation  to  the  roots,  espiecially 
as  it  increases  with  angle  of  attack,  is  discussed  et  length  in  Peference  1 
It  is  shown  that  it  is  a  good  approximation  for  lew  angles  of  attack,  I.e.,  3°, 
but  becomes  unacceptable  for  angles  of  attack  near  10^  especially  for  neutral 
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Etabllity, 

1.   The  General  Caee. 


zf 


Assuming  a  solution  of  the  form,  x  =  x  e    ,  vhere  z  ie  a  real  or  com- 
ples  constant,  the  Indicated  operations  on  the  variables  are  then  carried  out. 
The  stability  quartic  rfesulte  from  an  expansion  of  the  following  determinant 
that  follows  from  the  fact  that  the  equations  are  consistent  if  the  determ  nant 
of  their  coefficients  vanishes. 

Bear  longing  the  equations  by  variables, 
u 
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The  aseumptions  are  made  that  G^    sin  '^^'  f  is  negligible  since  near  horizontal 

flight  is  assumed,   and  that  aporoximately,   T     =  D   . 
'  o         o 

2  /  2 

The  determinant  of  the  coefficients   is  then,    UBin<3  k     =  S/i 

y       '   y 
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The  quartic  becomes, 

B.l  Z^  +  a  Z^   f- 

where, 


a^'      +  a,Z  +  a     =  0 
2  1  o 
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The  usual  methodE  for  the  solution  of  quartics  may  be  employed  to  8olve 
this  equation  as   it  stands.      The  conBtants  have  been  arranged  In  the  particu- 
lar faehion  shovn  for  reasons  to  be  mentioned  later. 
2.      The  Case  of  rietinct  Sets  of  Roots. 

Experience  has   shown  that  the  iisual  motion  in  flight  consiets  of  a  long 
period,    lov  frequency,    lightly  damped  motion,    called  the  phygoid  or  flight 
pat>^   oscillation,   and  a  short  period,   heavily  damped  motion,   called  the  rotary 
oscillation.     This  BUQgests  the  factoring  of  the  characterietic  eqiiation  Into 
tvo  quadratic  equations. 

The   len^ith  of  the  period  of  the  rotary  oscillation  has  been  found  to  be 
at  most  of  1-5  seconds  while  that  of  the  phygoid  usually  is  of  the  order  of 
30-60  seconds  for  most  conventional  planes.     This  eu^ests  that,   due  to  the 
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short  duration  of  the  rotary  oscillation,    the  velocity  is  eesentially  con- 
stant.    Thus  this  mode  is  governed  "b"  equations  A. 5  end  A. 6  vlth  the  terms   in 
u  eliainated.     The  solution  for  the  phygold  can  neither  ignore  the  change  in 
velocity  or  an^le  of  attack.     However,    cone id erst ion  of  equation  A. 6  ehova 
that,   due  to  the   If.rge  raat^nitude  of  the  factor y<^   diecussed  in  Appendix  A, 
an  approximation  may  be  made  that  the  terms  containing  this  item  are  large 
compared  to  the  other  terms.     Solving  this  equation  for    OC  ,    the  terms  in  (X 
in  equations  A.^  and  A.  5  are  then  ex  pressed  by  a  term  in  u.      The  phygoid  vrlll 
"be  solved  for,   as  regards  the  information  mentioned  previously,  by  considera- 
tion of  the  equations  enclosed  by  th.e  dashed   lines  in  the  equation  arrs.y  modi- 
fied by  the  approxirrEtion  diecuesed  above.     The  rotary  oscillotion  will  be  con- 
cerned by  the  equations  enclosed  by  the  dotted   lines   in  the  arxtiy. 
a.     The  Phugoid. 

Equation  A. 6  becomes,   nelgectlng  terms  not  involving yW  , 

hQ 
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Equations  k.k  and  A. 5   in  coefficient  form  are  now,   assuming  a  solution 


of  the  f  onn,   x  =  x  e  '        ,  ^ 
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which  gives  the  quadratic, 

3.2  Z^  +  a,Z  +  a  =  0 

1     o 

where 
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The  cuadi'atlc  has  the  roots, 
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Where  the  damping  factor,  E,,  is,  P>-,  =  -  ~  the  frequency  or  period  factor  iSj 


I,  = 
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1       */     o       14-1 
b.     The  Short  Period  OscillE'.tion. 

Nee^lectLng  the  velocity   terms,  -and  equation  A. If   in   its   entirety,    equations 
A. 5  and  A. 6  "become,    in  coefficient  form, 
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^liich  gives,   -ii-^   f  A.Z     -f  ys,     ,  =  0 
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Where,   A.   =    (J^        -  C      k" )   an'^,,   A     =    (j 

the  others  are,    '/'■  =  E^  n-il^  and  the  damvin^s  factor  is,  B.^  = ^.    the  frequency 
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'Before  any  further  evaluation  of  the  quantities  determined  up  to  this 
point,  the  quasi-Etatic  stability  criteria  are  developed. 
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C.      Qiiasi-Statlc  Stability 

Before  procoeding  vith  the  quasi -etc tic  stability  development,    it  may  be 
helpful  If  a  reiteration  be  made  of  the  specialisation  aeeumed   leading  to 
thle  developnent  compared  vith  methods  of  a  :nore  genei-al  nature  as  in  references 
1  and  k.     Both  of  these  reforencee  discuss  the  «pecial  cese  of  the  approximation 
to  the  roots  resulting  from  the  aseuraption  that  the  solution  consists  of  tvo 
pairc  of  complex  roots   suffici&ntly  separeted  in  nagnitucie  tc  wan-ant  a  factor- 
ing of  the  characteristics  quartlc     This  facile  solution  is  ueuil  in  most 
cases  but  may  not  bo  sound  for  unconventional  airciaft  or  even  conventional  . 
aircraft  where  additional  facte  must  be  considered.      These  are,    flexing  of  the 
component  parte  such  as  wln^^s  and   tail  surfaces  of  the  aircraft,   vhich  vould 
entail  additional  decrees  of  freedom,    large  changes   in  the  coefficients  of  the 
equations  due  to  rapid  depletion  of  a  large  fuel  eupplv  and  its  effect  on  the 
derivatives.     Many  of  these   latter  cone id orations  are  undergoing  exhaustive 
study  at  present  and  do  not   lend  themselves  to  any  concise  generalizations. 
Prom,    say,   a  project  engineer's  point  of  view,    especially  in  the   initial  stages 
of  design,    comcem  is  given  to  irformation  available  from  the  most  rapid  and 
inexpensive  data  available.      This   is  especially  true  at  present  since  specifi- 
cations,  the  military  in  pai-ticular,   require  employraer.t  of  the  finished   design 
t>irou£hout   large  ranjoe   in  altitu'^e,    speed  and  wei^jht.     The  cssumption  of  the 
presence  of  two  quite  widely  separated  tnodee,    however,    is  scvjid  for  most  conven- 
tiorial  aircraft  ajid  for  unccnventiCTial  aircraft  at  certain  phases  of  their 
flight  histoi'y.     The  oonditione  that  hold  for  the  phugoid  were  seen  to  be  a 
Blcnr  oscillation  of  long  period  and  wave   length  and  relatively  light  damping. 
The  rotary  oscillation  is  a  short,    heavily  damped  oscillation. 

In  regard  tc   etrmdard  methods,    the  cocfficlente  needed  to  determine  the 
CQnetar.ts  for  the  final  solutions  cf  the  equstions  of  motion  ane  obtained 
primarily  from  static  vind  tunnel  teete.     This  is  the  rule  fox'  the  values  of 


21 

equetiona  A.l  and  A. 2  ae  eeen  in  references  1,  3  and  k.     The  moment  coeffici- 
ents are  found  by  static  stability  coneideretiOTie  for  the  variation  of  the 
mOmer.t  witV  <X  .  See  reference  5"  A  dynamic  wind  tunnel  test  iB  used  to  find 
C   .  See  reference £  1,  k   and  6. 

q 

The  coefficiente  aB  they  hare  been  developed  in  thie  paper  vill  be  evalu- 
ated in  much  the  same  manner  as  regards  thoee  concerned  vlth  linear  forcee. 
Hovever,    the  method  of  evaluation  of  the   lo^rithmic  derivatives  vill  be  Il- 
lustrated.    As  for  the  moment  derivatives,    the  detemination  of  C       will  not  be 

m 

q 

discuEBed  aeeuming  that  it  may  be  found  by  the  usual  neihods.      The  static   longi- 
tudinal stability  criterion,   0         ,    is  the  quanitity  that  the  next  fev  pajce  wl. 
be  concerned  with.      This  q.uantity  Dr.    Scheubel  determines  in  a  quasi -static 
vay.      Qauci-£tatic   in  that  it  is  developed  frocn  assumptions  made  on  the  d-mamlc 
notion  of  the  aircraft.      The  tvo  concepts,    quasi-static  stability  at  equilibri- 
um and  at  constant  speed,    follov  from  the  fact  that  an  aircraft  has  tvo  dis- 
tinct phases  of  motion  following  a  disturbance. 

Cuasi-statlc   stability  at  equilibrium  of  forces   is  eeeentially  a  considera- 
tion of  the  equilibrium  ccnditionE  reached  at  a  tine   in  the  flight  history 
following;  a  distiurbance  at  which  the  nev  velocity  is  obtained.     This  new  veloci- 
ty is  the   initial  value  plus  the  incremental  increase  due  to  the  perturbation. 
The  cua£i-Btatic   stability  at  constant  speed  development  holds  fox^  the  rela- 
tively s>-ort  period  following  the  disturbance   in  vhich  the  velocity  is  assumed 
not  to  have  roaches   its  final  value. 

The  d:','Tiamic  response  to  an  elevator  deflection  is  disectod,    so  to  speak, 
to  give  those  two  quaci-Etatic  concepts.     The  aircraft  is   initially  in  steady, 
nearly  horizontal  straight  flight,      starting  from  this  steady  straight  flight, 
we  a£Eun3  that  the  elevrtor  angle,    ^    ,   bae  been  chango'd  suddenly  b-  a  certa  n 
amount,    d  ^  ,   and  we  ack  vh«.t  vill  happen. 
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In  the  first  inetant,  the  eleratcr  displacomejit  gives  e  chejige  in  longi- 
tudinal coraent  and  a  etnall  change  in  lift.  Thie  latter  change  ie  to   Enell 
it  is  neglected.   The  change  in  moment  dieturbing  the  equili"brium  gives  an 
angular  velocity  about  the  lateral  axis,  and,  due  to  this,  the  angle  of  at- 
tack is  changed  too.   A  change  in  the  angle  of  attack  gives  rise  to  a  change 
of  the  lift  coefficient  and  hy  it  a  change  in  lift.   Soothe  equilibrium  of 
forces  perpendicular  to  the  direction  of  flight  ie  dieturbed  too.  All  these 
effects  hold  during  the  initial  portion  of  the  period  folloving  the  distur- 
bance, i.e.,  during  the  first  few  seconds.   In  addition,  it  is  aseumed  that 
the  speed  remains  oEsentially  conetejit  during  this  period. 

This  then  is  the  citxjation  maintaining  for  quasi-static  stf.billty  at  con- 
stant speed  and  considers  the  change  in  force  perpendiicular  to  the  flight 
path  caused  by  the  initial  curvature  of  the  flight  path  and  the  change  in  mo- 
ments resulting  from  the  disturbance. 

Eventually  the  change  in  elevator  cnusing  the  change  in  angle  of  attack 
thereby  the  lift  coefficient  results  in  the  aircraft  attaining  a  new  steady 
state,  that  ie  a  new  point  of  equilibrium  of  forces  and  moments  at  a  certain 
speed,  V  *-   dV,  which  is  different  from  the  initial  one  by  dV.   Cuasi-etatic 
stability  at  equilibrium  concerns  itself  with  this  portion  of  the  flight  his- 
tory. From  experience  it  is  known  that  it  takes  an  aircraft  an  appreciable 
time,  normally  several  minutes,  for  the  speed  to  adjust  itself  to  a  changed 
elevator  displacement. 
1.  Quasi-static  Stability  at  Equilibrium. 

The  implication  of  an  initial 
steady  stoto  of  motion  necessitates  equi- 
librium of  forces  and  of  moments.  From 
Figure  k   this  is  seen  to  be,  "^ 
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L  -  W  cos  ai^=  0 
M  =  0 
Mr  »  0  . 


Having  started  fx-om  nearly  horizontal  flight  the  flight  path  will  remain 
nearly  horizontal  for  sufficiently  snail  changes  so  that  the  component  of  the 
weight  perpendicular  to  the  path  remains  almost  xmchanged.  From  the  second 
equation,  that  of  equilibrium  of  forces  perpendicular  to  the  path  of  flight; 

then, 

4(L  -  W  cos  if*'  )  =  dL  =  0,  since,  as  before  sin  Ir  =  iT  =  0. 
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where  the  subscript  E  indicates  an  equilibrium  consideration  with  the  aseump- 
tione  implies  above.  This  equation  is  neither  a  partial  nor  e  total  derivative 
from  a  tnathetnatical  point  of  view  and  caution  must  be  exercised  in  handling  it. 
If  the  new  state  of  motion  is  a  steady  one,  which  it  must  be, 

dC  =  M  -  M  =  0.  The  moment  coefficient  depends  on  0^  ,  V,  6    ,   and  the  angu- 

d0 
lar  velocity,  tt-  .  Since,  however,  it  is  assumed  that  the  flight  path  has  no 

curvature  at  the  new  steady  state,  and,  indeed,  in  actuality  such  would  be  the 

case,  the  angular  velocity  is  zero.  Thus  the  quantity,  q,  does  not  appear. 
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There  ie  then, 


ac 


or 


thUB, 


vhGje, 

C.2 


M  '  5^ 


dC.  = 


dcX 


he 


dS 


M 


d   In  V 


d   In  V  =   C 


M   ,     ~^M         d   In  Vn  _       ZZ}^  ^  f 


■'  =  -ST 


5} 


ST 

^^M       ^"M  d    In  V 

§^  ^  d   In  V         do< 


S=M         ^=M 


'E 


ho 
M 

dW 


--  X         --         d   In  V 
5^  ^  3    In  V         dO< 


This  then,    ie  a  new  etabllity  requirement  that  applies  to  that  phase  of 

the   longitudinal  moticn  deecrlhed  and  implied  from  the  previous  dlscusBion. 

This  requirement  certainly  is  sound.     For  an  aircraft  which  does  not  have  this 

qviasi-etatic   stability  vill  be  unstable   in  as  much  as  for  any  disturbed  state 

of  motion,    enf^ced  by  an  elevator  displacement,    will  show  the  tendency  to 

move   its  elevator  further  in  this  direction,    so  increasing  the  deviation 

from  its   initial  state. 

dCj^ 
,  Th^  usual  static   stability  criterion,   r-r-,    is  considered  either, 


^=H  ^=L 


dC^  6.0^ 

'   dC^  d«^  °^'  "  d^  ^^  "  ^'^*  '^^   quantity,  t-tt  ,  in  either  of  the  above  ie 

a  constant  throughout  the  usual  range  of  normal  flight,  i.e,  unstalled  flight. 
The  value  of  x  -  x.,  can  be  found  graphiclly  from  wind  tunnel  data  of  C^  vs. 
C-  as  explained  in  reference  5*  This  is  the  more  informative  quantity  since 
a  shift  in  the  center  of  gravity  is  more  readily  determined  throughout  the 
flight  hietoiT-. 

The  first  term  in  equation  G.2  is  very  similar  to  the  static  stability 
criterion  above.  However  for  comp8.rieon,  equation  C.2  is  written  in  its  en- 
tirety as. 


^J 


d<X 


(X  -  Xjj  Ij) 


IE     ^"^  'E 

where  x^j-r.  is  a  nev  neutral  point  that  takes  into  accomit  variation  of  the 
neutral  point  vith  velocity,  i.e.,  compressibility  effects.  First,  It  is 
noted  that  variation  of  lift  coefficient  with  angle  of  attack  differs  "between 
the  two.  This  is  due  to  the  fact  that  the  ordinary  determination  of  the  change 
of  lift  coefficient  with  angle  of  attack  is  accomplished  in  wind  tunnel  tests 
where  the  stream  velocity  is  kept  constant.   The  results  are  then  baaed  on  a 
constant  dynamic  preeeure,  wViereas  the  present  development  is  not.  The  dif- 
ference between  the  two  is  about  I-30/0  and  is  sumed  negligible.  The  quantity 
x-,|_  moves  with  respect  to  velocity,  dynamic  pressure,  deflection  of  the  fuse- 
lage, twiet  of  the  wings  and  horizontal  stabilizer.  A  comparison  with  the 
movement  of  the  usual  static  stability  neutral  point  will  be  made  only  vith 
respect  to  variation  with  velocity  and  dynamic  pressure.  Figure  5  shows  the 
variation  with  velocity  of  three  different  aspects  of  stability  at  sea  level. 
These  are:   (a)  static  stability  only,  (b)  quasi-static  stability  assuming  the 

zero  lift  moment  coefficient,  C^.  ,  to  be  -0.02,  and  (c)  quasi-static  stabili- 
se 
ty  with  zero  lift  coefficient      assumed  zero.  Figure  6  shows  the  comparisai 

of  the  same  quantities  at  altitude.  Here  the  decrease  in  the  stability  margin 
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for  qvasi-etatic  etability  is  much  leee  duo  to  the  presence  of  the  density  in 


the  dynanic  preseure  in  the  denominator  of  the  addlticnal  term,  i.e.,  C_  in 

d  In  V 

E 


d  W 


2.  Cuaei-Static  Stability  at  Constant  Speed. 

Eecalllng  the  dlBcuseion  of  the  condltiOTis  holding  following  the  imposing 
of  a  small  elevator  deflection  on  the  steady  state,  it  vrs  seen  that  in  the  flrE 
small  interval  following  the  change  the  speed  does  not  change  appreciably.  Thu« 
the  quasi-static  etability  at  equilibrium  concept  does  not  hold  for  this  phase 
of  the  disturbed  motion.  It  was  further  seen  that  the  equilibrium  of  forces 
perpendicular  to  the  direction  of  flight  is  disturbed  too.  The  equation  of  mo- 
tion for  this  direction  shove  what  happens. 

If  a  /arc^   is  exerted  on  a  body  moving  in  a  certain  direction,  its  path 
vill  be  curved.   See  Fig\are  7-   If  an 
aircraft^ the  ijotion  may  become  rather 
complicated,  for,  due  to  the  curvatxire 
of  the  f  lifght  path  in  a  vertical  plane, 
the  component  force  of  the  weight  per- 
pendicular to  the  path  changes  too. 
Eeetrictior.  to  email  disturbances  close 
to  horizcaital  flight  allows  neglecting 

the  changes  in  the  weight  component.       "  Figure  ?• 

For  instance,  a  change  in  ajigle  of  flight  path  of  15  gives  only  3  c/°  change. 
The  centripetal  ac:;eleration  is  V  rf^  •  The  equation  of  motion  is  therefore, 


and 


^V  11=   d(L-  Wooer) 


dL 


dL 


■[: 


Sc, 


>xd0x  dt 
^^dt^  - 
he. 


^f^ 


-  'r7*«  -s?^^>prs 
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therefore,  ^^  ^^y  o 

dS       ,^^L  ,  «      *'^L  ^   X  o  ir     s 


m 


and,   multiplying  by  — ^  and  noting  that  —  =  f^ — ~ ,   there  is 


or, 


C.3 


/4dq(l  -^  )   =5^dCX 


d£ 

do( 


The  term,   — -j ,   may  be  neglected  since  5—  is  small  compared  to  5-—- 

and,  also,  the  ueually  Ipr^e  magnitude  of  aa  makes  it  negligible  compared  to 
unity,  /his  is  also  the  comncn  approximatiosi  as  seen  by  the  neglecting  of  X 
and  Z     in  references   1,    3  QJid  k» 

A  I  T  1 

Thus  there  is  the  relE.tion,  r^    =  vrr-  -rj-  ,  resulting  from  a  considera- 
tion of  the  forces  perpendicular  to  that  of  flight  in  the  first  instant  follow- 
ing a  distux'bance. 

This  curvature  of  the  flight  path  or  angular  velocity  has  an  influence  on 
the  equilibrium  of  moment  about  the  lateral  elxIs.   So  a  consideration  of  the 
equilibrium  of  moments  along  this  curved  flight  path  involves  the  change  In 
elevator  deflection,  the  change  in  angle  of  attack  and  this  angular  velocity. 
The  effect  of  velocity  is  omitted  because  of  the  assumptions  of  this  phase  of 


the  motion.   There  is  then. 


Sw.      d 


t«. 


or, 


he 


M  ^   M  dq 
SoT  '  5q~  do(  y 

This  expression  is  similar  to  the  one  found  for  the  equilibrium  of  a  steady 
state  of  flight  enforces  by  a  certain  elevator  displacement.   It  contains  the 
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etalDility  terra, 

PL  5l  ^^M      ^^     i2 

^•^  d«     y  ■  S^      ST    do^ 

called  "quael-ctatlc  stability  at  conetent  speed".     A  comparison  of  this  quan- 
tity vlth  the  usual  etatic  stability  is  seen  folloving  a  further  evaluation. 
The  angular  velocity  causes  the  angle  of  attack  of  the  horizontal  tail 
plain,  P(  ,,    to  change  by  the  a.mount, 


d0 
h  =  ^  °"^'  dt  ==   ^^  ;V2 

therefore, 


,  o*     dt  ^h         ,  .       d0   ,^  V 


^^h    \ 


dq    5T2   ■  ■ 

This  change  in  tail  angle  of  attack  gives  riee  to  an  increase  in  lift  coef* 

ficient  of  the  horizontal  tail  plane,  C_  ,  and  by  it  a  decrease  in  moment  about 

^1 


tbe  lateral  ajcis, 

1 

^=H=-\>I 

and,    having  noted 

its  dependence  on  kngular  velocity,   there  is, 

dq 

d0(^        ^3/2        dq              dcX^      h   ^    0^ 

This  quantity,  Cj.  ,  may  be  determined  by  a  dymanic  wind  tunnel  test  de- 
q 
scribed  in  reference  1,  or  by  use  of  the  curved-flow  teclinique  referred  to  In 

reference  6. 

hC^    d   I 
The  term,  r-—  t^  I   ,  added  to  the  static  stability  is  alvays  negntive, 

since  C   is  a  damping  derivative.  Thus  this  is  an  increase  in  static  stabili- 

q 

ty  and  the  neutral  point  corresponding  to  this  quaei-Btatic  stability  is  be- 


hind t>^e  —  "--  .---.-   -.-^-..-.^     .      ...   .       K 


hind  t>'e  one  for  ttatic   stability  by  about,    —  ,    vhere,   K  =   C        J_ 

This  apparent   improvement  over  wind  tunnel  static  stability  ie  seen  to  be 
larger  for^an  aircr-aft  of  small  vlng  loading  at  low  altitude  than  for  one  of 
the  high  wing  loading  at  high  altitude. 


The  addlticnal  term  added  to  the  ueual  etatio  etabilitT  neutral  point 
loolcE  £ue?lciously  like  the  ctick-fioced  maneuver  point,   eepecially  because  of 
the  etatement  of  the  situation  and  the  aeeumptione  ntide.      This  ie  indeed  the 
caee  and  ccd.  be  verified  by  conpariecn  vlth  the  mcjieuver  point  found  in  ference 

u. 

i 

Concerning  the  concepts  of  quasi -stati^:  stability,    it   is  seen  that  both 
are  improvements  over  vind  tunnel  stttic   st^ibility.      Hovever,    three  facts  should 
be  borne   in  raind  ccfnceming  their  effect  ."nd  their  use.     First,   both  cont^iin 
the  uEual  etc  tic  stability  plus  an  additional  term       Secondly,    their  applicatior 
and  use  only  hold  foi-  specific  conditions  perculisr  to  different  phases  of  the 
flight  hietor-r       Lastly,    for  a  specific  aircrr.ft  they  vary  differently  with  like'- 
parametorc,    i.e.,    density,   but  all  parameters  do  not  effect  both 

Dr.   Scheubel  pointed  that  in  Germany  during  the  lo.Bt  var  certain  of  their 
fighters  shoved  a  marked   improvemont  in  longitudinal  stability  over  that  pre- 
dicted by  vind  tunnel  tests  at  medium  altutudes.      tince  these  were  of  the  kyj- 
5CO  raph  class,    this  vas  explainable  by  the  above  quasi-static   stabilit;"-  consid- 
erations.     He  said  fuithur,    however,    thct  pilots  coniplB.ined  of  poor  stability 

f 

at  high  Epoece  at  high  altutudes.   In  as  much  as  both  concepts  still  predict  en    ■ 

improvement,  greater  for  equilibrium  and  less  for  constant  speed,  the  aeronau- 
tical engineers  vere  greatly  puzzled.  Further  study  acooiaited  for  this  by  rea- 
scn  of  flexing  of  the  wings,  the  tail  plane  and  the  fuselage.   The  r.c   ountablli- 
ty  of  these  phenomena  in  stability  studies  ie  very  complex.  For  the  longitudi- 
nal case  a  good  approximation  made  be  obtained  from  grouping  their  overall 
effect  as  a  shifting  of  the  neutral  point.  The  seriouenesB  of  their  effect  can 
be  seen  from  mention  of  two  simples  aspects.  First,  upon  the  sudden  applicetior. 
of  an  elevator  deflection  or  a  change  in  speed  as  in  a  gust,  the  ©L5.sticity  of 
present  designs  at  high  speeds  causes  a  twisting  of  the  vlng.  This  causes  the 
angle  of  attack  to  very  epanvise  resulting  in  a  shifting  pf   the  center "of 
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f 

pressure.  The  same  applies  to  the  horizontal  tr.il  surface.   Secondly;  e. 
shifting  of  nev-tral  point  vertically  f.-om  flexing  of  the  fuselage  coupled  vith 
the  very  large  drag  forces  at  high  speeds  produces  an  apprecia"ble  deeta'bilizing 
moGient. 

That  the  problem  ie  complex  analytically  ie  easily  seen  from  merely  a 
superficial  study  of  the  present  day  literature  as,  for  instanfce,  reference  7- 
Dr.  Scheubel  pointed  out  the  experimental  difficulties  in  a  description  of  an 
attempt  in  Germany  to  determine  the  deflscticn  of  a  horizontal  stabilizer  vith  , 
speed.  This  involved  a  complex  optical  system  employing  an  arrangement  of 
mirrors  reflecting  a  slit  light  source  the  movement  of  vhich  was  recorded  by  a 
motion  picture  camera.   In  addition  to  this  the  test  aircraft  was  required  to 
dive  along  a  path  deterninod  by  an  ar.-ay  of  floodlights  to  insure  perfect  lin- 
earity of  flight  path." 
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I),   ivaluetlon  of  the  rerivatlvee. 

Sone  of  the  quantities  to  'be  evaluated  are  iienticsl  to  those  appearing 
in  standard  referenctc.  M-thode  for  the  determination  of  theee  will  "be  men- 
tioned but  not  discuEsed  in  detail. 

Since  all  quantities  appear  in  the  general  solution,  the  items  vlll  be 
diecuseed  in  the  order  thr.t  they  appear  in  the  coefficients  on  page  11. 


'lA 


^M. 


£t?tic  wind  tunnel  teste. 

lyn^.mic  wind  tunnel  test  by  the  whirling  arm  or  oecllL=i.tion  me- 


thods mentioned  in  reference  1  or  by  the  ciirved-flov  techjiique  referred  to  in 
reference  6. 


^  In  D 
d  In  V 


Then, 


A   ^lot  of  C-^  vs.  V  can  be  obtained  from  static  wind  tunnel  tests 


D  =  C^  P  ~  S 


=  ==f 


dV 


T 

2 


1 

V 


2) 


V  dD 
D  dV 


d  In  r   ^  ^  S 


d  In  V  "  ^  In  V 


Thufc  this  varies   f -  om  a  vr-ilue   of  2 


at   zero  M  uo  to  about  3   f-r  M  =    1.0. 

^    In  C^ 

Theee  aje  con- 


^    Ir.  T  ^    ,  "    "•  ^T 


sidcred  t0c-;cthsr  because  they  are  found 
in  nuwh  the  same  manner.   From  wind  tun- 
nel tecte  C^  vs.  V  are  foun-i.   Ingine  Fi£;Txre  8 
Tianufactui-urs  usually  submit  data  of  T  vs.  V  for  specific  in  stf- list  ions.  Those 
derivatives  are  then  determined  as  in  Fi^gures  9  and  IC.   It  is  pointed  cut  that 
Figure  10  is  plotted  for  a  constant  angle  of  attack. 
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Figure   10. 


Vo^m 


g— r — ?r  .   A  plot  of  Z     v^,   V  can  be  found  from  wind  tunnel  tests.   The 
derivative  is  then  found  as  shovn  in  F'gure  11.    , 
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nrv  =  ^  ^M 


FisvLT-e  11. 
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E.      Evaluation  of  the  Longitudlnel  Modes  Using  GuBsl-Static  Concepts. 

1.     The  Phugoid.  • 

The  damping  factor,  K, ,  depends  on  the  coefficient  a  of  eqtiation  B.2  and 
this  quantity  is, 


a  In  D   S  In  T 


®1  "  "D  ^d  In  V  "  B  In  V 


At  relatively  lev  velocity  the  approximate  values  of  the  derivatives  are, 
<  ?^  y  =  2,  ?  j^  y  =  -1  (for  a  .let)  and  ^  ^  y  =  0.  Thus  the  damping,  is, 


^1 "  "2  Ts^^A  =  -  2  s  ^^'^^  ^^  ^^s^^- 


At  higher  speeds  the  effect  of  the   last  quantity  in  E^  will  l3e   felt.      This 
vill  decrease  the  damping  since  it  will  overcome  the  increase  due  to  the  other 
quantities. 

The   period   factor,    I,,    depends  on  "both   the  coefficients   of  equation' 3. 2. 
The  coefficient,   a  ,    using  the  quaai-etatic  stalDility  at  equi  libriiira,    equation 
C.2,  "becomes, 


*«M 

'-L^"       S    Ir.  V   ' 

c„ 

d« 

iL 


The  period  factor  nov  become^s, 


..; 


o  ^    In  G 


^M 


li 


F 


'MjK 


1  /^    In  D       o    In  Tv2 
IT  ^S    In  V   "  S    In  V^ 


=  c^V^^ 


d    In  G 
d    In  V 


,S    In  D       S    In  Tx2 


F       1     ^o    In  V       6    In  V 
■  IT  d    In  Ct   _ 

'2^    nnr> 


Again,   at   low  speeds   the  approximations  c^n  be  made  thnt, 


i^r<<  1 


>I 


and, 


o( 


=    1 


"K 


^^ 


So  that, 
Kow, 


d    in  C^ 
^1=   ^lJ^-  -"     h    In  V 
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/    a  In  G_  '   r 

S 

In  C^  ' 

T     p   /  o  1        J:t            1  '' 

u 

V 

Ln  V 

The 

period  is  then, 

^    2n       2it 

V 

s 

and 

the  wave' length  is, 

.   L  =  V.T  =  ^ 

P       P    /    d  ]^ 

c/ 

/-  1  . 

ii 

V"  ^  1X1 

V 

2. 

The  Short  Period  OscilLation. 

8 


The  damping  factor  is,  E^  =   -  ^  (C^       -  C„     k  ) . 
Approximate  values  for  these  quantities  are, 

^L       =  5.C,        C^     =    -0.7,  and,  C^^ .   k^  =    -p   . 

The  damping  factor  is  the,  P.  =  ■4'  5  -  (-5)  =  "5  which  indicates  heavy 
damping. 

,  The  period  or  frequency  factor  was, 


/  ^ 
^2  W'  ^ 


Since,  -  ^  (C^^  -  C^^  k^)^  =  .  ^  (:^^)2  ,  ^  c^^c^^  k^  -  ^  (C^^  k^)^ 


2  2  2 


yU  =  -  C.,   k  ii  -  G.-  C-   k  ,  the  eum  is. 


-  IT  ^^Lol^^  ■  ^  ^L„,  Si  ^Y  "  ^  ^^M  ^v^^  ■  ^M  ,  ^l/*-      ^^SroMVixiS   this  quan- 
tity, 

-  C„   k^  M  -  r  (Ct   +  c,,  k^)^  . 

X 

And,  using  the  values  given  "before,  it  ie  seen  that  the  item  in  the  parenthesis 

is  zero.  At  least  the  value  in  the  parenthesis  is  very  small  and  it  "being  to 
the  second  pover,  divided  by  four,  in  addition  to  the  presence  of  the  Large 
f&otor, yii   ,    in  the  first  terra  allows  neglecting  it  in  comparison  with  t>-e  firs 
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is  now,  Ig  =  V^Lc^  ^"^  "  \^  K/^  '    ^°  ^^^^' 


term.   Thue  there  xe  u<^w,  -"-^  -  ^  v^-r  _,  v^   ^/  "-^r/M 


CaJ  = 


I 2 


^cV2  JL 


Co     /c     ' 


1 


The  period  is  the^i,  T^  =  -^  =  2«  / -^^ r  ^ut,^  =  r?»  ^"^ 


So  then  the  period  "becomes, 


wc— li^y ■  ■  ■  ■«—«* 

.2 


S    y2 


Frotn  the  foregoing  it  is  seen  that    depends  on: 

1 


yf 


l(     center  of  gravity  location: 
2)     Altitude: 


^u*y — r 

3(     Wing  loading:       -^{^^i      • 

U)  Size:   -VW  .        ^-,,-1 

/i" 

5)  Maes  Concentration:   y  -J 

6)  Velocity:  |  . 

The  decrement  ie   seen  to  he, 


2«  Eg       (=L^  -  °M^,  '^y'^ 


2       V=L^  (='  -  ^'  I'y/. 


Although  thie  analysis  is  specialized  and  the  modes  of  motion  are  of  such 


frequency  and  damping  that  they  do  not  preeont  much  ccncom  in  design,  the  ap- 
plication of  the  quasi-ote.tic  concepts  has  "been  ehown.   In  regai-d  to  the  im- 
portemce  of  the  mod.es,  reference  8  points  out  that  the  importance  of  the  phu- 
goid  should  not  "be  minimized  too  much. 
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The  evaluation  of  the  coeff Icicnte  of  the  non-factored  qufirtic,    equation 
B.l,    page    11,    in  terras   of  the  qinei-Etatic   ete/blllty  concepts  vill  not  be  cr.r- 
ried  out.     The  coefficients  have  been  arranged  In  n   form  for  easy  eubetitution 
of  theee  quantities  and  a  glance  vill  Bhov  by  vhat  degree  they  can  he  simpli- 
fied by  such  a  substitution. 

It  is   impoi'tant,   hovever,    that  a  criterion  for  the  separation  of  the  quar- 
tlc  Into  modes  of  slow  and  fast  frequency  be  eetabliehed.      If  the  ratio  of  the 
imaginajry   pai-t  of  the  roots  describing  the  phujjoid  to  that  of  the  short  period 
oscillation  is  snail  conpawred  to  unity,    then  the  approximation  is  acceptable. 
Fron  ths  simplified  frequency  factors  used  in  the  previous  analyeie,   '^.nd^    it 
is   pointed    out,    the  terms  retf.ined  are  the  meet   important  ones,    it  c-?n  be   seen 
that  this  criterion  is   primarily  dependent  upon  the  fhctors  z  -  x^.,  and   ii    in 
the  nxmerator  of  I^.      This  follcvs  since  the  other  terras,   althcui:h  the;"  tnay  var^ 
are  rels.tively  constsjit  vhile  the  distance  bet^jeen  the  center  of  £pr?.vity  and 
the  neutral  point  may  be  anywhere  -from  O.GO5  to  O.C5.      The   lar^je  variation  of 
the  density  factor  is   sho^nn  in  /ypendix  A. 

The  presentatiCTi  thus  far  has  dealt  with  the   ttick-fixed  case   onb'.      The 
degree  of  freedom  about  the  elevptoi   hin^e   line  is  considered  next. 
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F.   Introduction  cf  the  L levator  Force. 

In  Section  A   it  ie  recalled  that  the  equ,?tion  of  longitu'^.inal  motion  due 
to  the  degree  of  freodoto  cbcut  the  elevator  -wsa  locked,  i.e.,  it  Implied  that 
a  change  of  elevator  angle  ^^B  Impreeeel  and  holl  faet.  Kov  this  i*eetrictlon  ie 
removed  and  a  result  obtained  in  the  form  of  a  modifier tl on  to  the  original 
eolutioT:  of  the  ehcrt  period  oecillntior..   This  Ie  the  etick-free  cr.ee.   It  ie 
interest in<3  to  note  here  that  the  most  cotrnicn  texts  on  this  euh.iect  treat  this 
case  quits  differently,  reference  h   rules  cut  '^■ny   modification  of  the  phugold 
mode  and  soIvcl  the  reniainint^  equations,  dropping  terms  in  u  from  -.n  assumption 
of  constant  speed,  foi-  tvo  short  perioc  -nodes.  One  of  the  Letter  is  dieregardec 
as  the  flapping  cf  the  elevator  about  its  hinge  line,  while  the  other  ie  diecue- 
sed  extensively  as  the  mode  causin^^  the  "porpoisinG"  motion.   Peference  1 
jfrites  off  the  freeing  of  elevators  as  a  modification  of  the  modes  by  e  reduc- 
tion  of  the  tail  effect  by  reason  of  the  ving  downvaeh  and  a  possible  contribu- 
tion to  the  velocity  moment  derivative,  K^,  by  I'eaeon  of  some  constant  moment 
due  to  elevator  mass  or  a  spring  attachment  or  friction. 

The  equation  of  motion  resulting  from  this  degree  of  freedom  consiate  of 
equating  the  moments  resulting  f^'om  freeing  the  elevator  to  the  inert ial  tern, 
or  moment  of  inertia  of  the  elevator  times  the  angulax  ac  .elei-ation. 


/ 


Figure  12, 


Thus, 

F.  1  rn  1  — '^-^ :r — ^  =  Mr 

dt^  S 
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From  the  figure  it  is  seen  that  a  dovn  (positive)   elevator  deflection 
causes  a  nose  do^<m  (negative)   -noment  about  the  aircraft  center  of  gravity.      It 
was  neglected  for  the  stick-fixed  case  but  will  nov  be  included.     This  quantity 
a^ears  in  equation  A. 6.     There  are  effects   in  directions  coneietent  with  the 
other  equations  also  but  these  are  assuraed  small. 

It  is  also  reasonable  to  assixne  that  there  ie  an  increase  in  moment  about 

the  hin.-^re   line  of  the  elevator  due  to  the  change  in  an^le  of  attack  of  the  tail, 

Ae  is  pointed  out  in  i*eference  h,    this  ehoulcf  be  reduced  to  as  near  poesibl© 

ae  possible.     A  typical  fighter  of  the   last  World  War  had  a  value  of  C         25  o/c 

tbat  of  C,   ^    .      C,        may  be  made  to  be  smr^ll  by  protier  design  and  many  methods 
h^  h<x       •'  ^    y     ^  ^ 

of  reducing  it  are  discussed  in  references  h  and  9-      The  incremental  angle  of 
attack  of  the  tail  itself  is  small  for  the  portion  of  the  flight  histcr;/  con- 
eidered  here,   and  even  that  is  considerably  reduced  bv  ving  downwash  depending 
on  the  horizontal  tail's  position.     The  C         term  vrlll  bo  neglected  here,  but 
certainly  should  be  included  if  design  prevents   its  being  made  small       Any 
large  value  of  C,        will  detract  from  the  pilot's   "feel"   of  the  aircraft. 
Derivatives  of  the  hinge  moment  with  respect  to  rates  of  change  of     O^     anri  5* 
are  nelgected  also. 

Ae  in  reference  h,   any  modification  of  the  phu^^oid  from  f-'eetng  the  ele- 
vator will  be  acsuraed  to  be  of  no  consequence.     Thus  the   improvement  to  the 
etick-fixed  approximation  will  only  consider  the  short  period  osclllat  on,, 
which,   a£  mentioned  before,    assumes  constant  velocity.      Thus  this  equation  will 
neglect  all  changes  with  respect  to  velocity. 

The  second  term  in  KU   is  the  coriollis  force  which  results  from  reference 
of  the  angular  motion  about  the  aircraft  center  of  gravity  to  the  center  of  gra 
vitv  of  the  elevator.     The   Last  tern  is  the  mass  moment  effect  due  to  the  angu- 
lar acclcration  about  the  aircraft  center  of  gravity.      So  F.l  becomes,    letting 
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*=L.v2 


A  =   -  ^-^  P  h-  £.    C     and  non-d imenBlone.il zing, 


a     i 
e     e 


2       ^   ^  _      ^  „  -^^.,.2 


e-S       .    r  .  "^e    Po^  .       d  S-  "^e^fe  ^h       ^e^     d^(  OC       ^)     ^ 


A  <     d(f-)2         "  ^  S     *  d(^)-^                A  t;             d(f  )^           '      '. 

Kov  lat,                         P  p       ■->                       '          p      l/^ 

ml  m     i     27^                         m     r  S"^" 

_     e     e  e     e                                  e     e 
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•"-    •    .2 


e 


i"-  i~ 

e  e 


then  F.l  "becomes, 

F.i        t  (1  :  ^)ir  -  t  bxa r  +  t  (1  +  ^)6<    t  *s    5*  =  o. 

B  e    y^  e  e  " 

Thi£   eqxiation  along  vith  equations  A. 5  end  A. 6   of  Section  ;!^   gives  tVe   evs- 

tem  which   ehove  the  modee  of  motion  for  the  stick-free  ccee.      The  complete  s-'e- 
ten  of  equations,    includintS  A.U,    should  he  coneldered  for  a  rigorous  analvsie 
of  the  modee  of  motion.       IIowGver,   as  vas  diecussed  pr-oviouely,   the  short  peri- 
od oscillation  presumes  conetajit  velocity  and  this   l£  ago  in  assumed.      Thus  there 
is  the  following  set  of  equations. 


A. 5  r  -  ^T  <^  0  =  C 

A. 6  f    -  ra    JT'  OC    -  n   «    -  ra  .  iA.O<  -  ™gM  T  =  ^ 

F.l     t.d^)**--  t  3^f  t_(l'-^)0<  t^S    -S  =  Co 

222 

'      q         My'  0(  M^      7   '    .  "S         Mg     y 

Aeeuming  a  solution  of  the  foi-m,   x  =   x     e^  ,    there  is  the  detenalnf  nt, 


ko 


z     -  n^z  -  m^^ 

0 

t  (1  +  ^)z^ 
© 

=   0 


or, 


2  '^  2  I 

(t     z     +    1)   (z-*  -  a^  z     -  m.,        lAz) 


where  recs.ll, 


which  is,    regrouping, 
„3 


z2(C^^-^^)    ^    .(-m^l^)/*    '    t^[-^0 


z^   +    z2{C, 


-     (1    -^)2''(C, 


(X        <1 

z)mr/U-Btim<-uC^      z  >.  =  0 


m  )    -  zm, 


Lm] 


ijyu-  a^<,B5^v-j_ 


4" 


or,    eliminating  the  root,    z  =  G,   and  letting,    f,(z)  =    z     +    z(C^ 

i  i-«0(        <i 


-  tn 


-  V  Iv/^ 


=  c 


)j     3iithere  ie. 


and,    f  (z)  =   (1  .  ^)(z^  +  C^      z) 

F(z,t^)   =   f^(z)        t^j^.2   ^^^^j    ^   mg^f^Cz)]  =   C. 

If  now  the  case  is  considered  for  the  raaee  density,  m   ,    ij^rge,  ^^.ich  Is 

usually  the  case  except  for  glidera,  then  the  qus.ntitv,  t  ,  essentially  vaniehe. 

£o,  approximately  then,  f,(z)  =  0.   The  entire  expression,  F(z,t  ),  may  be 

X  e     '^ 

analyzed  from  the  point  of  view  of  using  the  roots,  z.,  of  f,(i.)  In  the  evalua- 
tion. 

The  effect  of  freeing  the  elevators  then  will  be  the  detertnlnotion  of  the 
damping  factor  and  frequency  of  the  short  period  oscillation  from  a  factor,  ^Iaz, 
to  be  determined. 

Beviewing  the  situation,  the  characteristic  eque.tion  obtained  from  the  de- 
terminant gives  a  fifth  degree  equation  in  z.  Elimination  of  the  null  root 
gives  F(z,t  )  which  is  fourth  degree  in  z.  I^eference  U  then  solves  this  equ-^- 
tion  for  two  paire  of  roots  which  give  the  two  short  period  oscilLe.tions  .-e 
pre viouElv  nientionad.     - 
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There  is  thon,    in  mathematical  termB; 

t     =  0  giving,  f,(z)  =  0  where, 

e  -I- 


z  =    z. 


then 


Sf 


dt. 


sr       ^sy^^c^'i^ 


e 
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f'fv-      ) 


Kow, 


but, 


GO, 


and. 


(l^^)z^,.    (1       ^)r^^,   -S^3^ 
=    -^g/* ^ 


^1  =^  ^1   '    ^^1 


p         p         o 
.2   =  p-    .   q  -:-   2i  E^  I, 


^^Lc^-  \^   +     ^/'^^cIvA   ^'l 


"l=     ■ 


■^; — r 


2  2 


(c,   -  n  y 


(C       -  m^)^ 


(xly/*  T 


=    -m 


c<'v/* 


therefore. 


,2       ^2 


dz 
dt 


-   -^jf/^ 


(1      ^)(F.^  -  I^  .  2i  B^  I^    -    (B^  -    iI^)C^^)    -  C^J^ 

21;      ^  2il      -  ?K 
1  -^1       -   1 

(1   !   ^){2i  B,    I,    ^    II,   C^      )  ♦ 


=     -^Sy^ 


1      1  1     L 


(X 


2i   I. 


-   i 


L.  1  -^ 


21   I. 


=    -m 


s> 


(1  -  u)(2R,   -  0,      ) 


but, 


1  „^) _ > 


k2 


2E^    ■    C.     =    -C^      +  rn      -  0^=01 

2^1  ^  ^  ^L  -  -2-  -  'l^\  -^  i  -  t-  ^  — JC" 

=   -  ^     (C       -  m  ) 
2  L^        q 

th  crefore,  ^p 


(1  '  a)  [m^  l^y.^  !^  (^  -  m^)]-   C^^3^ 


2^1 


but  now,     /AZ  =  rr—  t       .     Coneidering  the  term,    -m^  jut       , 
e  ^ 


_  ,2   ^1/2  V.2 

-C        /i  m     i     S  ^      k 

'f/*  "e  "  -Cv,  r  u  m  S.    C. 


but, 


2      .,.  ^"^L.        ..2 


V"   ^3/2  .    .  ^  p  1 


2     W 


o 

thuE,  ^C^ 

\         ,2 
\  r     m     i     r. 
m     .   4-  o  ^        e     e     h 

"™5/*^  =  -5c: 

n     ^     c^     r. 

coneequer^^y,  ^  " 


(i^:^)[^/.    ,^(=,^.V]'  V^ 


It  ie  seen  that,  since  both  C    and  C   are  negative,  the  real  part,  which 

S  \ 

indicatee  damping,  ie  negative  and  thiiiB  this  motion  is  damped  for  this  accep- 
tance of  sign.  However,  the  value  of  C^  might  go  positive.  For  inotance,  re- 

8 
ference  9  ehowe  that  a  horizontal  tail  design  using  an  KACA  66-009  airfoil  eecti  ii 


^3 
with  flap  0.30  chords  in  length,   an  air  gp.p  of  C.OO5  chorda  end  with  the  hinge 


lir.e  at  a  positon  of  C.50  flap  chords,   reducee  the  C         to  ^ero  between 
0(  =«   -8°  to  0<    =    '5°,   while  G^        is  positive.        This  reference  also  ehowB  that 
C         nay  30  positive  depending  on  the  use  or  non-uBe  of  abov*»  oal?-nce.     The 
use  of  a  horn  hale-nce,    in  general,  appears  to  mr.ke  C         go  slightly  positive 
alao,    the  oagnltiKle  "bein^  ve-^y  Gmall. 

Thus  this  mode  of  motion  ae.j  be  positively  damped  or  possess  moderate  di- 
vergence.    The  fi^equency  is  given  by. 


e  e 


where  the  teiin  ne^jlccted  is,   (j-   -  m  )  ,  and  is  considered  small  ne   in  the 

previout  diecuESion.  The  terms  in  the  brr.ckete  are  negative,  the  last  term 

positive.  The  total  value  of  the  quantity  is  not  too  large. 

It  is  difficult  to  etrike  s.  concrete  ccoparison  of  the  mode  of  motion 

solved  for  here  and  thoee  found  in  reference  h.      This  solution  does  not  solve 

the  equation  by  a  f©ctori2ation  but  by  the  aesumption  of  the  quantity,  t  , 

being  large.  In  general,  the  roots  solved  for  by  Dr  Scheubel  are  not  t'-e  s-^me 

as  thoee  causing  the  porpoising  motion  described  in  reference  h,    since  there  it 

is  indicated  that  the  nature  of  the  porpoising  mode  is  a  function  of  C.    -^nd 

C.   .  Eere,  it  is  dependent  upon  C,  -  ,  -rrr     and  C   .   It  is  not  to  be  ira- 
h  ^  h  £    '     6  b  ra 

plies  that  this  is  not  such  a  ty^pe  of  moti'On,  rather  that  it  is  not  the  same 
implicitly. 

Dr.  Scheubel  neglected  the  quantity  Z         ,  it  ie  believed,  due  to  the  fvTct 
that  in  Germany  It  wos  the  practice  to  design  the  controls  such  that  it  was  re- 
duced to  a  m  nimujn.  He  discussed  theT  h^jidling  qu/ilities  of  ?  very  large  four- 
engined  flying  boe.t  in  v/hich  the  stick  forces  per  g  vers  very  emr  11  without 
power  boost,  the  exact  values  of  which  are  not  remembered.   In  regard  to 


kk 

I 

lateral  control,   however,    it   Is  remeniberod  that  he  etated  thct  one  pereon 

could  naintain  rtreight  and   level  flight  with  "both  engines  dend  or  one  side 

with  (Xie  foot  on  the  rudder  ccntrol,   a^^in  vithout  power  "booeted  controls. 

Inclucion  of  tlie  C.        torm  would  pL-ce  it  In  the  solution  of  f,(z)   eo 

that   it  would  detract  or  add  to  both  the  damping  and  frequency  terms  as  it  wne 

negative  or  poeitive  in  the  rc.tio  of  r; Furthennore,    in  os  much  aa  down- 

"h^J 

wash  would  multiply  the  value  of  the  vs-riable,    0(    ,    which  rffects   it  "37  a   fac- 
tor less  than  unity,    this  would  fui*ther  decreiiSe  the  ratio  itself. 
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G.     Effect  of  Elcvp.tor  Impulse  Type  Forcing  Functicti. 

F.ei  tea's  tine  J    it  is  recalled  thet  the  stick-fixed  cr  ee  involved  the   impor- 
tance of  en  elevetor  dieplBcemor.t  or  o   ^^uat  and  the  eneuint;  motion  w^e  nnnly- 
zed   in  the   li::ht  of  the  elevrtor  being  locked  in  poeition.      The  etlck-free  caea 
alloved  the  elGv:tv-r  to  flr^et  freely  a"bo\;.t  its  hinge   line.      Furthermore,    the 
equations  of  notion  were  aseumed  hono{.;onQoue.      An  additionri  1  cree   ie  now  con- 
sidered  in  vhich  the  cquetiona   ir.volve  a  forcing  function  i^nd  further  this 
forcing  function  is  an  elevrtor  impulse,    n.n  elevr-tor  deflection  of  short  dum- 
tion.     This  cace   iE  analogous  to  the  incorporation  in  the  eystem  of  :'n  on-off 
type  autopilot. 

An  elevator  impulse  as  shown  in  Figure  13   is  imposed.      This   impulse  ie 
given  the   Gymhcl,    I  g    ,    and  ie, 


r  •  t^  =  /^  s  dt 


^  The  first  phase  following  the  impulse, 

thtt  of  the  short-period  oscillation,    is 
COT.eidered  and  thus  there  are  the  eqiiatione,  Figitre  13- 

The  effect  of  the   inpulee  on  the  oquf.tion  of  motion  perpendicular  tc  the 
wind  axis   ie  negligihle  for  conventional  .^ircrr?.ft  elnce  the  only  force  crerted 
ie  a   lift  force,   which,    due  to  the  short  duration  of  the   impulse,   does  not  act 
long  eroui^h  to  boaffective.      This   ie  not  so  fo     talllees  c.ircrr.ft,    however, 
since  for  euch  aircraft  the  deflection  of  the  elevrtor,    vhich   is  rictually  a 
portion  of  the  main  lifing  foil  trailing  edge,    changes  the  wing  angle  of  att-=ck 
and   caueer  nn  appreciable  effect   in  this  direction. 


K6 


The  effect  of  the  impulee  on  the  equation  of  moraente  c&ueee  on  Initial 
angular  velocity  which  will  be  deeicnated  oo^ ,  but  no  anguL'^-r  diepL?, cement. 
If  there  le  no  dieturbance  in  flight  path  angle,  which  is  nearly  correct  for 
conyentional  aircraft,  this  initial  angular  velocity  will  coneist  eeeontially 
of  solely  a  time  rate  of  change  of  angle  of  attack.  This  le  eo  becruee  of  the 

mass  dietribution  of  a  conventional  aircr?.ft  compared  tc,  spy,  r<   tAilleee  air- 

^  -  s   ^ 

craft.  Thus  Initially,  at  ^  =  0;   >»'  =  y*  =  0;   «  =  C;  o(  =  cOq 

Solving  A. 5  for  ^  ,  substituting  in  A. 6,  there  is  then  for  the  ehort  per- 
iod oscillation. 

The  solution  for  0(    ie, 

^  1  ^/         ^^  o* 

(X  =  e   (A  COB  I^   B^  sin  I,  ^   ). 


Imposing  the  first  initial  condition  gives,  A  =  0.  Differentiating  o(  then, 

(p(   =  UJ^=   e  ^     ^^1  'l  ^^^  ^1^   ^1  ^1  *^°^  ^i'^^  • 

The   second  initial  condition  gives,     B,  =  -r=—       .      The  Principle  of  Impulse  and 

■^1 
Momentum  states  that  the  change  in  angular  momentum  about  any  axis  equ:\le  the 


anj_,ular   imt>uleee  of  the  applied  forces  about  the  axis,    or, 

y 


mi     Lo     =   nr  S{t)      eo  that, 
y       o  6 


,.  Lie;  Srf.,/N, 

"^  0  =  ^  ;2 ^,  ^^^" 

mi  C 
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/>• 


V2..2 


C         TT   S  t 
■^TT-    i       - 

<^3/^      y 


f 


__c       ~    -^   P^i 


-5      i^       y^         -       '     V2        -^      '    ''^/'"     '^ 


y 

This  evalxjatoE  the   constant,    3   ,    so  thrt  consequently. 
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<X.^    I 

1 


3         Bin  I-   ^ 


From  equation  A. 5  there  ie,    X^  -  C^    0(     =  c.       — = —  I-     e         ein  I>.     Inte- 
ger ^t<       ^1       ^  ^ 


gratiiHj, 

(G.i)  r= 


B^'f 


■S        2 


-^  (I^^  ein  I^^   -  I^  coe  I^i:  )   ^    ^ 


Since  at,  "T^^  0,     Jf*'^:  C  there  is  finally, 

G.l     r=  -^ —^     I.       1  -  e   '-       (cos  1,^  -  r" 


ein  l^"t    ) 


This  initial  phase  of  the  response  to  an  impulee  tvpe  forcing  function  is 
depicted  qualitatively  in  Figtire  Ik.      In  equation  G.l  the  vcluea  of  B.   and  I. 
are  the  same  ae  thoee  found  in  Section 
F.,  page  2,6.     Thus  the  transient  is 
heavily  damped  and  typical  values  show 
it  tc  have  reached  the  cteady  state  in 
one -ha If  periods. 

Examining  further  the  steady  state 
of  G.  1, 


Figure  Ih. 


^   = 


'^/^"l, 


A'u's 


n  ■  n 


and   lettin^::,   G        L.  =    -C_      (x  -  x,.!,.)   there   is,      y    =    ■    '•■       i  '  If  the  de- 

nominator  is  snail,    the  rec  pcanse  is  very  sensitive  to  the   impulse.      Dr     Scheu- 

bcl  cited  an  exajnple  of  neutral  stability  vhere  the  steady  Btat«  then  becomes, 

C  I^ 

^  ''  since,  x  «=  x.. 


-  C       "^ 
\ 


!1_   y«A 


Kow,   C 


-r^       ^/-.        ana,    -  C       —  =  '  ^ 
og     ,,3/2  n 


V  =   U50  mph;   D  «  270   ft.'^;   I,.  =    -1°;    t     =   0.2 


^h  ^h        1 
q/*      ^^t       S^       y^ 


—  .      Using  values  of, 


sec;   u  =    IOC,    the  response  is, 


1*8 


So 


-/P      .V 


*^.  =  -  -^  •  V^  •  -T75  •  ^J  •  ^a  -  •'■•'   ^  ^  ^  1^  '^  -  "-"-^S  =  5.7° 

which  is  appreciable  for  euch  n  small  impulee.   It  is  pointed  out  that  this 

value  of  clerator  cffcctivGnese,  —rzz  >    correepcTi:  s  to  rn  elevator  area  to  to- 

'      a.  ^ 

tal  horizontal  tail  area  ratio  of  C.6  according  to  an  omplrlcc.lly  determined 
curve  in  reference  k. 

The  second  phase  of  tho  motion  following  the  impulse  ie  thnt  of  the  flight 
path  oscillation.   Thuc  eqiir.tions  A.k   and  A. 5  are  ueed.   These  should  be  modi- 
fied ae  before  for  the  variable,  O^  .  However,  it  will  be  o.eeumed  that  the  de- 
rivativo,  5— r — rr  ,  is  approximrtely  zero  and  will  be  neglected. 

Fron  the  pact  analysis  of  the  two  modes,  the  reL  tive  magnitudes  of  their 
periods,  ar.d  the  najiitudes  of  their  damping  factors,  it  can  "be  i-issumed.  that 
the  response  to  the  inpultc  of  the  rotary  motion  has  reached  the  Etia.dy  state 
condition  by  the  time  the  flight  path  response  takes  effect.   So  the  Initicl 
conditions  are,  at,  'X  .^  C :  S'  -    if   ,    and  u  =  0,  a,»n  d  J**—  o. 

Again  the  forcing  function  creates  nc  accountable  effect  either  along  the 
wind,  axit  or  perpendicular  to  it.  Thus  a  ;:ood  app2-oxim;'tion  is  to  vork  vith 
homogeneous  eqiiationt.  So, 

A   solution  for  <r  ie,  ^=  e  "   (A^  cos  l^Z   >-   B^  ein  I  'ZT  ).  The  first  ini- 
tial  condition  gives,  at  ^  =  0,  A^  =  y*  .   The  second,  B_  =  -  7^  ^   since, 


y  =  C  =  B  A   -  I,j  B^.   So  nov. 


h9 


0.?. 


K"^ 


r  =  dr^  e  -       (coc  I,^  i'    -  ~  Bin  I^  C-  ) 
The  inclination  of  the  flight  peth,    S*"  ,   rar.y  "be  cefinod,    using  Fi^-ura  1^, 


ae, 


vu       1     ^.h 

^  =  V     dt 


i,:i^l±... 


3 


r,, 


ano,    puttinG  this  in  terns  of  th- 


corimon  time,    ^     ,    tincc,   t     =M — rr-  -=  -~     vh 
,^    l-i : '      a     y         V / 


ere. 


/.  =,^  ^'^' 


eo  thr.  t, 


Subttltutin-  for   ^  and   inte3rat:ng, 


^£0  -^^ 


.2        .,2 


?r.  I. 


2         2 


..S 


"2  '  I--^        li'  \  J.  -  o         coe  I  -J  6 


■2     ^    K^£ 


The  maximum  altitude  gained  foUoving  an  Imyulee   is  then,   eeeuraing  T     and 
o    In  C,  •  2 


h    Ixi  V       ^^11^ 


h  = 


I2    ^» 


^. 


^'=-^^ 


o    in 


n     00 


f? 


gT- 


ThiE  rceult  rne.y  be  obtained  by  a  GliGhtly  different  Qpproc.ch.'     Coneic^erine  equr 
tio-n  A,  5  obove, 


u  = 


o    In    J 


=l(2   ^  5ir-^)       '2  *   ^2 
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Differontiatinc  G.2, 


r  =   -  y..  e  '-         ^  ^     -   em  I,  L   . 


00  1 

2 


Therefore, 

"r. 

e  '"     sin  I  ,'6    . 


Goneidering  maximum  u, 


J 


and,  neglecting  ener^  lo£S  due  to  dampins, 

H  -^  h  •  —  =  constant 

V    \  '  V         'V  V 

keo  that,      d(h  h-  —)  =  0  or,   dh=-—     -^=-— ^u     fejid,   finally, 


U 


II  V  I         I  V^  ^=0 

I  Ah  I  =  "7"  1^1   "=  T"         •   ^       which  it  ae  before. 

Usin^  the   same  fi^jui'es  ae  "before,    the  mt;jcinnan  t^Jiin   in  altiti^e  ver-  ap- 
proximately "becomes, 

_h  =   '^  -  ■  '     :<  0.0995  =   95c  feet 


C 


(H^ 
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DlscuE8iC3n  and  JonclurioniB 

The  quasi -static   etcbllitv  concepts  based  on  equilibrium  of  forces  snd  tno- 
raents  at  varicue   phases  of  the  fli^^ht  hiEtor^r  of  cjn  aircrift  following  r:  dis- 
turbance appeal'  to  bo  eound  for  an  accept-inae  of  the  deacription  of  the  flight 
hietory  statsd  earlier.     NQglectlng  addltionrl  effecte  euoh  -e  eV^eticit^'-,    these 
Qonctptc  account  for  the  api)£rent  incresee  in  etr.biljtv  duo  to  compress ibi lit v 
effects  associated  vith  th(-   aircraft  drgi3,    lift  an'   roomer t       If  the  v  lues  of 
these  quaei-Etotic   Etability  criteria:   coula  be  rscert'  Ined  Individu-^lly  and   In 
their  entirety,    oreeumably  they  vculo    include  intex-ference  efx*ects.      Aa  it  ie, 
no  American  literatui-o  known  to  the  writer  discusses   lonjitudln-'l  str-bility 
comparisons  zt  low  and  hi^h  speeds       Br.    Scheubel  discuse^e  the  coraxv^rieon  with 
the  requite  cho\fl^.  in  F  ^J^iros  5  ^^^    6.      Thia  ^*as  aponrentl:'-  b-ised  cc.  results  of 
flight  tests  or  the   inclusion  of  the  additional  factors   In  the  quasi-etatl'c 
criteria  analytically  determined.      This   L  tter  would  necessarily  neglect  in- 
terference effects.  ^ 

For  the  special  case  of  sepc.mtlon  of  the  motion   into  the  phu^old  ."aid  shor- 
period  modes  good  a^jreement  v/lth  etrindard  developmentc  was  had  neglecting  the 
additional  factors  mentioned  above. 

The  modification  of  the  phu^oli  mode  due  to  freeing  the  elevators  \ros  ne- 
glected as  Ie  usually  done.      The  reculting  equrtion  wps  analyzed  for  the  spe- 
cial cace  of  the  aircraft  nr^as  density  factor  of  1  r^e  m.-agnitude.     This  assum- 
ption vae  used  to  approximate   only  one  pe.ir  of  the  complect  con'u^te  roots 
Comparison  with  standard  references  shoved  this  motion  to  be  approiimatelv  the 
"porpo*Siny"  modi  althcu^;h  not  implicitlv. 

Of  some   interest  is  the  GermDn  use  of  the  quantity  S"^,    square  root  of  the 
wins  area,   ae  the  non-d imensi ana li zing  length  term  in  convereion  of  forces  and 


53 

momente  to  coeff iciente.  ActualXv,  Dr.  Scheu-oel  did  not  point  thlB  out  as 
having  wide  Gertnpn  uea^e,  "but  rfithcr  ae  "being  U8©d  nt  rarTn8t.''.dt  Technical  In- 
Btitute.  It  m.^ht  be  pointed  out  thet  the  concept  of  r.epect  rrtio  vould  be 
l08t  by  ure  of  thle  len^h.  Thie  raicht  be  better  than  the  nebuloue  definition 
of  taean  aerodynanic  chord  present  in  current  Amcricrn  aerod^mamJc  litem.tvire 
For  inetance,  reference  U-,    euppoaedly  a  fairly  CiCcredited  text,  defines  mesMn 
aerodynamic  chord  as  1)  an  average  chord,  for  a  finite  ving  of  vi-.ryino  section, 
vhich,  vhen  multiplied  by  the  average  section  moment  coefficient/  dynamic  pres- 
sure and  wins  area,  gives  the  noment  for  t^e  entire  wing;  2)  the  chord  of  an 
imaginary  airfoil  which  throughout  the  normal  flight  iT^nge  has  the  ecme  force 
vectors  as  the  three -dimensional  wing  and  ceji  be  determined  from  vind  tunnel 
teste  or  a  gi'ayhical  analysis  which  assumes  c   special  shape  to  the  ving  pL'^in, 
and  3)  a  chord  on  which  all  th<?  forces  and  moments  acting  on  the  surface  crji 
be  represented  as:  acting  and  on  which  the  pitching  moment  coefficient  is  in- 
variant with  lift  cocff ici^int.  Dr.  Scheubel  pointed  out  thr.t  the  speciol  rea- 
son for  the  ucc  of  this  quantity  was  the  fact  that  the  ratio  of  the  squrre  root 
of  the  wing  area  to  the  tail  length  was  found  to  ri.nge  between  0.8  and  125  ^^nd 
the  average  value  found  for  some  ninety  aircraft  v/ae  1.01.   The  rdvcnt^^.ge  of 
approximating  this  ratio  to  \mity  resulted  in  the  adoption  of  thie  length. 

The  extension  of  such  a  method  to  the  lateral  or  as-rmmetric  case  ie  ([le- 
cussed  in  gone3ral  terms  in  Appendix  B. 

The  development  for  the  longitudinal  case  as  a  whole  does  not  consider  all 
the  factors  bearing  on  the  problem  by  any  means.   It  was  pointed  out  by  Dr 
Scheubel  that  this  was  a  good  approximation  from  the  prelect  engineer's  view- 
point and  must  necessarily  be  improved  as  the  design  develops.  l''or  inat-ince. 
the  effects  of  downwash  and  the  influence  of  the  poverplont  have  a  L-^.rge  ef- 
fect on  the  damping  term.   In  any  case  it  is  a  slightly  different  view  of  the 
problem  and  its  voir  greets  on  ttet. 
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•*  APPENPIX  A 

The  taaes  niim^er,    aa   ,   as  mentioned  in  reference  1  vr.e  firet  eu^ jested  "by 
Glauert  for  uec  in  the  nan-limonQianalizinc  of  the  equatione  of  motion       Thot 
it  £leo  follovE  from  diraeneione-l  an£\ljDie,    e.^^.,   the  ?i  Theorem,    ie  well  knovn 
Tue  to  ite  extenrive  ucg   in  the   litoreture,   a  diecueeion  of  it  here  ie  cxil'"  ws-r- 
rantod  "by  the  sli^^htly  diffor-ent  definition  of  it  in  this   psper  ojnd  the  fact 
that  Gonslderatione  of  its  magnitude  fox-  uae  in  the  approximatione  in  thepaper 
neceeeitate  recalling  ^rtiat  factors  bear  on  ite  magnitude. 

The  difference  in  definition  liee  in  the  use  of  the  mean  vin^  chord  ae 
the  represontrtivc  length  in  references  1  and  h.     Furthermore,   the  raoee  or  den- 
sity numher  ueed   in  this  paper  is  tv/ice  that  used  in  theee  references 

The   quantity  written  aa, 

from  vhich  it   is  seen  to  depend  on  altitude,    O    ,    on  vin^  lending,   W/S,   and 
onxeight,    l/W.      The  tendencies  with  these  parameters  v.nryins  are  as  follows; 

a.  Altitude   -  As  altitude  increases,  y*^  increases.      An  altitude  increfise  from 
sea  level  to  i<-0,000  feet  causes  an  increase  inyt^  of  from  1  to  U 

b.  Win^^  Loasffing  -  Since  it  is  to  the  3/2  power,   an  increase  of  ving  loading 
of  1:2  gives  an  increase  iniA.  of  1:2.83. 

c.  Weight   -  Ac  weight  increases,   x*    decreases  as  (]/W)^.      in  this  case  a   1:2 
increase  in  veight  give  a  1:0,7  decrease  in  Ai   .      Typical  vclues  ofyU.    ere: 

Lerge  gliders  at   low  altitude:    v-   =   5-IC 
Small  fighters  at  high  altitude:   y*    =   5OC 
large  coraaercial  plrjiec:    lA  =  5O-IOC 
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APPniDIX  B 

LiscusBicn  of  the  Extension  of  c.ucsi-Strtic  St;  "bill ty  to  the  Lateral 

or  Aey.mnfctPlc  Motion 

The  possibility  of  tho  oxtensicn  of  the  quasi-etatic   stability  concepts 
to  the  asyrnraetric  motion  of  an  aircraft  seems  to  "be  remote  insofar  r>s  the 
sane   line  of  reasoning  be  folloved.      It  will  be  recalled  that   in  both  quasi - 
staiic  Etabilit^/  dotertainations  consideration  vra.s  made  of  forces  in  a  single 
linear  direction  and  of  taomonta  about  a  point  in  one  plane.      In  effect  the 
ratio  of  chanses  in  pararaeters  affoctinQ  both  eyetema  lead  to  the  quasi -et<a- 
tic  stability  concepts.     Of  fundamental  importance  vac  the  fact  that  consider- 
ation of  those  changes  did  net  affect  the  additional  linear  force  equs-tlon. 
£uch  a  correspondence  does  not  hold  in  tho  lateral  case. . 

/.asuraing  a  conventional  aircraft  it  has  been  found  that  the  modes  of  mo- 
tion in  the  lateral  caee  are  four  in  nurr^ber.     Following  the  same  line  of  rea- 
soning as  vac  done   in  tho  loncitudinal  case;   a  disturbance  could,   all  factors 
considered,    cause  a   linear  displacement  laterally,    i.e.,   along  the  y-axis,    a 
rolling  displacement;   and  a  yaving  displacement.     References  h,    IC  and   11  shov 
that  the  solution  of  the  lateral  equations  results  in  a  characteristic   quar- 
tic  which  gives  roots     \  ^     \  ^,      /^,,   and     X  j^  .     In  solving  for  these  roots 
all  nerjlect  the  presence  of  products  of  inertia,    that  is,   they  aseume  that  the 
principal  longitudinal  axis  of  the  aircraft  in  is   line  with  the  flight  path, 
and  assume  level  flight  or  tan  0=0.      The  moment  compoaed  of  three  components 

given  by: 

n  h'^  ^1^  ^1^ 

I  ^^=  B^  c-"     ;   p^=  P^e     ■"     ;     r^=  E^e     ■" 


I.  -  -.' 
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III    ^.  =  3^e     -^        Bj^e  ;     p^  =  P^e     -^     -^   ?^e  ;  r^  =  I^e^         :  j,  e 

vt.Gre  the  quantitiee     ^  ,    p,   ;  nd  r  r:re  in-'icctive  of  sideelip,    ••nctilrr  veloci- 
ty of  roll,   end  rn^jular  velocity  of  yew  roepectively.       The  coefficients  F   , 
P.   and  3     r.re  dependent  on  initial  conditione,    onV  four  being  e.rliitr&ry,   the 
remainder  are  dettrmined  fiom  these.      These   indivldur.l  mctiaae  vrLll  "be  con- 
sidered separately. 

Considering;  the  first  group,    the  root,    A   ,  "is  small  rjid  ueur  lly  positive 
giving  light  divergence.     This  is  the  so-called  spiral  motion  and  either  p.  emrl 
dieturbance  i/hlch  ic  ^ui-c  sidoelio,   roll,    or  yr.w  will  develop  an  interdepen- 
dence of  both  rolling  and  y^ving  rnomer.te  on  "both  sideellp  angle  and  yawing 
velocity.     Following  quasi-Gtatic  reasoning,    sinue  the  period  ie   long,   this 
mctior   indicatee  the  eetabliehment  of  a  unique  reLrtionehip  between  initinl 
dieturbance  incremer.t£  and  resulting  forces  and/ or  momentB  in  equlllbriian 
after  a  long  period  of  time.       This  does  not  follow  since  either  aileron,   rud- 
der or  side  clip,   due  to  a  gust,    sav,   nv?.y  cauee  the  motion  and  all  fncete  of 
latera,l  aotion  are  involved  in  the  resulting  motion. 

The  Ecconcl  group  ie  chsracterized  by  a   large  root,     A  ,^   usually  negr.tive. 
B^  and  F     are  sme'.ll  in  comparieon  vith  P     so  t'':.?t  the  motion  is  heavil-"  damped 
and  ie  oecentially  pure  roll.     To  knovai  reference  cor.eidere  this  motion  of  ony 
importance.       Peference   L2  which  includes  the  effect  of  finite  product   of  In- 
ei-tia  and   initial  inclination  of  the  flight  prth  excludes  "i^  mention  of  it  oIeo. 
Eeferencc   13,   which  ic  eji  inveetigation  of  unertlefactory  lateral  st-  bility 
OTi  an  actual  current  cor.bat  aircraft,    does  net  mention  it.        This  evidence 
coupled  with  the  feet  that  the  motion  involves   coaly  one  degree  of  freedom 
ce*tAinly  dooc  not  lend  itself  to  any  quasi -et;- tic  stability  criterion. 

The  third,  group  is  concemed  with  the  all-important  oscillator--  motion. 
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For  fixed  controlc  thie  coneiEte  of  "both  a  rolling  motion  and  a  yawing  motion. 
It  tec  crji  "oc   induced  "by  S'riy  one  or  nny  com'binction  of  the  posBTole   lateral 
di&tui'bancee. 

Jrotn  the  foregoing  it  ic  ccncludod  thr.tvTiique  representations  in  the  quaei' 
eta  tic  senbe  cannot  oe  deduced  for  the  aBjrmjietric  motion.      It  appears  from  the 
conclueionE  of  refcrenco  13  that  a  poc^i'bilit:'-  m  ght  be  a  c^sideration  of  the 
decree  of  frs^om  alaout  the  rudder  hin^e  line.      In  the  investigation  outlined 
in  this  report,   rudder-fized  considerations  were  ruled  out  in  that  acttxal 
flight  tests  ehoveu  good  ae^reonent  vith  predicted  oscillations  up  to  M-ich  num- 
bers corree ponding  to  a  opeed  of  h^t  ciph  (density  altitude  at  this  speed  vas 
not  specified).      Insofar  a&  a  gradual  1  deterioration  in  damping  with  speed 
resulting  in  rudder-fxeo  tc-cts,   a  variation  of  the  rudder  hinge  momentB, 
both  0,        and  C,  ,   vith  Mach  number  is  inferred.      It  is  also  stated  in 

this  repcr-t  that  this  trend  has  been  observed  in  high  Mach  number  tests  of 
other  contraol  surfaces.     Amplification  on  these  latter  raentioned  tests  is 
not  ua"de.     '2\to  facts  nust  be  pointed  out  here,_  howeypr.     Firsts   the  above  men- 
tioned possibility  is  not  in  tho  quasi-static  sense  of  Dr.   Scheubel's  approach. 
Second ly,   reference  13  based  its  computed  values  at  high  speed  on  an  extra- 
polation of  a  determination  at  an  airspeed  of  210  miles  per  hpixr  assuming  a 
ccjnstaiit  value  of  the  stability  derivatives  throughout  the  speed  range 

A  fur-ther  complication  in  the  establishment  of  quasi -static  criteria  fcr 
the   lateral  case  are  the  inclusioa  or  exclusion  of  tho  effects  of  products 
of  inertia.     These  effects  are  dic-cuesed  in  detail  in  reference   1^+  and  men- 
tioned in  reference   12  for  a  specific  aircraft  conf i(_:ur'-tion.      The  effects 
of  their  inclusion  was  considered  negligible  in  reference  13,   but  this  cor- 
clueion  must  be  made  weighed  in  the  li^ht  of  the  remarks  madie  above. 
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